E2-2000-217 



Very High Multiplicity Hadron 

Processes 

J.Manjavidze a f\ and A.Sissakian b 
(To appear in Physics Reports) 



"Joint Institute for Nuclear Research, Djelepov Lab. of Nucl. Problems, 141980 Dubna, Russia, 

E-mail: joseph@nusun.jinr.ru 
Tel.: (09621) 6 35 17, Fax: (09621) 6 66 66 

6 Joint Institute for Nuclear Research, Bogolyubov Lab. of Theor. Physics, 141980 Dubna, Russia, 

E-mail: sisakian@jinr.ru 
Tel. (09621) 6 22 68, Fax: (09621) 6 58 91 



Permanent address: Institute of Physics, Tbilisi, Georgia 



1 



Abstract 



The paper contains a description of a first attempt to understand the ex- 
tremely inelastic high energy hadron collisions, when the multiplicity of produced 
hadrons considerably exceeds its mean value. Problems with existing model pre- 
dictions are discussed. The real-time finite-temperature 5-matrix theory is built 
to have a possibility to find model-free predictions. This allows to include the 
statistical effects into consideration and build the phenomenology. The questions 
to experiment are formulated at the very end of the paper. 
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1 Introduction 



The intuitive feeling that hadron matter should be maximally perturbed in the high energy extremely 
inelastic collisions was a main reason of our effort to consider such processes. It was a hope to observe 
new dynamical phenomena, or new degrees of freedom, unattainable in other ordinary hadron reactions. 
This paper presents a first attempt to describe the particularity of considered processes, to give a review 
of existing models prediction and, at the end, we will offer the field-theoretical formalism for hadron 
inelastic processes. 

Thus, considering hadrons mean multiplicity n(s) as a natural scale of the produced hadrons 
multiplicity n at given CM energy y^s, we would assume that 

n » n(s). (1.1) 

At the same time we wish to have 

n << n m ax = y/s/m, (1.2) 

where m ~ 0.1 Gev is the characteristic hadron mass. The last restriction is introduced to weaken the 
unphysical constraints from the finite, for given s, phase space volume. We should assume therefore 
that s is high enough. 

The multiple production is the process of colliding particles where kinetic energy is dissipated into 
the mass of produced particles 0. Then one may validate that the entropy S accedes its maximum 



in the domain (1.1) since the multiplicity n characterizes the rate of stochastization, i.e. the level of 
incident energy dissipation over existing (free) degrees of freedom. 

There is also another quantitative definition of our reactions. Let e ma x be the energy of the fastest 
particle in the given frame and let E be the total incident energy in the same frame. Then the difference 
(E — e max ) is the energy spent on production of the less energetic particles. It is useful to consider the 
inelasticity coefficient 

E Srnax -, £max ~ /-i o\ 

K = = 1 < 1. (1-3 

EE~ v ' 

It defines the portion of spent energy. Therefore, we wish to consider processes with 

1-k«1. (1.4) 

So, the produced particles have comparatively small energies. 

This property may be used for experimental triggering of our processes. Indeed, using the energy 
conservation law, 

ra(l-/c)>l. (1.5) 



Following (1.2) we will assume that 

m 

(1.6) 

Therefore, the kinetic energy of produced particles in our processes can not be arbitrarily small. 

Using thermodynamical terminology, we wish investigate the production and properties of compar- 
atively 'cold' multi-hadron (mostly of 7r-mesons) state. We would like to note from the very beginning 
that we have only a qualitative scenario of such states which may be produced and the review of the 
corresponding why's may be considered as the main purpose of this paper. Although at the end of 
this paper (Appendix K) we will describe principal features of possible field-theoretical solution of our 
problem. 

It should be noted the absence of any authentic experimental information concerning discussed 
processes. Moreover, actually the hadron inelastic interactions with a set peculiar to them of unsolved 
theoretical problems will be considered. Nevertheless we suggest to work in this field in spite of these 
difficulties because the system with extremal properties may be more transparent since the asymptotics 
always simplify a picture. We would demonstrate this idea and will try to put it in the basis of developed 
theoretical methods. 

Absence of experimental information about so high inelastic hadron processes is the consequence 
of the smallness of corresponding cross sections. Besides this it was unclear for what purpose the 
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experimental efforts should be done. We would like to convince the reader that the discussed problem 
is interesting and important. For instance, we will discuss a possibility that asymptotics over n may 
replace in definite sense the asymptotics over y/s. A short address to experimentalists will be given in 
concluding Sec. 4. 1.2. 

We hope that the paper would be useful both for theorists and experimentalists. By this reason 
the main text of this paper will contain only the qualitative discussion of the problems and results. 
The quantitative proof, formulation of pure theoretical methods etc. are added in the Appendices. 
Considered extremal problem is a good theoretical laboratory and is described in the Appendices and 
theoretical methods may be applied for other physical problems. 

We would like to point out that the special technique was built for the problem discussed, see Sec. 2. 

- Having the very high multiplicity (VHM) state it is natural to use the thermodynamical methods. 
We will offer for this purpose the real-time S- matrix interpretation of thermodynamics. It can be 
shown in what quantitative conditions it will coincide with simpler canonical imaginary-time Matsubara 
formalism. We will give also the generalization of the real-time finite-temperature perturbation theory 
in the case of local temperature T = 1/(3 distribution, when (3 = j3(x,t). This will allow to use the 
thermodynamical description if the system is far from equilibrium. 

The particle spectrum in the VHM region is soft. It is just a situation when the collective 
phenomena should be important. To describe this phenomena, the decomposition on correlators will 
be adopted. The origin of this decomposition lie in the Mayer's 'group decomposition'. In multiple 
production physics this decomposition is known also as the 'multi-component description' It based 
on the idea that the multiple production process may include various mechanisms. 

In Sec. 3 we will investigate model predictions for VHM region. We would like to note two main 
conclusions: 

- Existing multiperipheral-type models are unable to describe the VHM region. 

- The infrared region of the pQCD becomes important even if the constraint (1.2) is taken into 
account. 



2 Qualitative inside of the problem 

In Sec. 2.1 we will try to formulate the phenomenology of our problem, i.e. the way the VHM processes 
may be described and what type of phenomena one may expect. The importance of thermodynamical 
methods will became evident and we will offer in Sec. 2. 2 a general description of the corresponding 
formalism. 

It is important to note that we may classify the possible asymptotics over n. We will find that 
there exist only three classes of asymptotics. This will simplify consideration definitely restricting the 
possibilities. 

In Sec. 2. 3 we will use the thermodynamical language to give a physical interpretation of these 
classes of asymptotics. We will see in result that in our choice of the VHM final state this should 
lead to reorganization of multiple production dynamics: we will get out of the habitual multiperipheral 
picture in the VHM domain. 

Moreover, one may assume that the semiclassical approximation becomes exact in the VHM do- 
main. This naturally leads to the idea to search such a scheme of calculation which depends on the 
choice of final state. Quantitative description of this idea may be realized as is described in Appendix 
K. 



2.1 Phenomenology of VHM processes 

The VHM production phenomena includes two sub-problems. First of all it is the dynamical problem of 
incident energy degradation into the secondary particle energies and the second one is the description 
of the final state. 

We will start discussion in Sec. 2. 1.1 from the second part of the problem to explain that the 
statistical methods are essential for us. 
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In Sec. 2. 1.2 we will try to outline at least qualitatively the main mechanisms of hadron production. 
The peculiarity of hadron production phenomena consists in the presence of hidden constraints, the 
consequence of local non-Abelian gauge symmetry. The constraints may prevent thermalization and 
the incident energy dissipation is confused in this case. Just the 'confusing' effect is dominant in the 
hadron multiple production processes if n ~ n(s). 

The expected change of dominant mechanism of hadron production is discussed Sec. 2. 1.3. It is 
important that, in spite of hidden constraints, the system may freely evolve to define the VHM state. 
Such a VHM state should be in equilibrium. Formal definition of the 'equilibrium' notion will be given 
in Sec.2.2.2. 

The problem described contains small parameters (n(s)/n) << 1 and (1 — k) << 1. To have the 
possibility of estimation of contributions in accordance with these parameters one should include them 
into the formalism. This becomes possible if and only if the integral quantities are calculated. So, the 
multiplicity n is an index only if the multiple production amplitudes a n (j>i, ■■■,Pn) are considered. But 
the cross section a n (s) is a nontrivial function of n. We will calculated by this logic mostly integrals 
of |a„| 2 , living out the consideration the amplitudes, see also Sec. 2. 2, where the first realization of 
this idea is offered (the naive attempt to realize this idea one may find in ||). This is a general 
methodological feature of our consideration. 

2.1.1 Formulation of the problem 

The multiple production cross section o~ n (s) falls down rapidly in the discussed very high multiplicity 
(VHM) domain ([O]) and for this reason the multiplicities n ~ n max are not accessible experimentally. 
At the LHC energy fi(s) ~ 100 is valid and we will assume that n ~ n(s) 2 ~ 10 000 is just the discussed 
VHM region (n max ~ 100 000 at the LHC energy). We will explain later why 

n~n(s) 2 (2.1) 

is chosen for the definition of VHM region. 

Generally speaking, having the state of a large number of particles, it is reasonable to depart 
from an exact definition of the number n of created particles, their individual energies Ei, momenta 
etc. since they cannot be defined exactly by experiment. Indeed, for instance, full reconstruction 
of kinematics is a practically impossible task because of neutral particles, neutrinos, the more so as 
n ~ 10 000 is considered. We suppose that nothing will happen if n is measured with An ^ accuracy 
since (An/n) << 1 is easily attainable in the VHM region. Besides it is practically impossible to deal 
with theory which operates by the N — 'in — 4 (~ 30 000!) variables. 

Artificial reduction of the set of the necessary variables may lead to a temporary success only. 
Indeed, the last thirty years of multiple production physics development was based on the inclusive 
approach H, when the measured quantities (cross sections) depend on a few dynamical parameters 
only. But later on the experiment and its fractal analyses shows that the situation is not so simple, 
also as, for instance, for the classical turbulence. So, the event-by-event experimental data shows that 
the particle density fluctuation is unexpectedly large Jq| and the fractal dimension Df is not equal to 
zero 

We know that if the fractal dimension is non-trivial, then the system is extremely 'non-regular' 0. 
So, Df ~0.3 for the perimeter of Great Britain and Df ~0.5 for Norway. The discrepancy marks the 
fact that the shore of Norway is much more broken then of Great Britain ||. It is noticeable that the 
fractal dimension Df crucially depends on the type of reaction, incident energy and so on. 

It is evident that one may choose from N = 'in — 4 an arbitrary finite set of variables to characterize 
the multiple production process. But the fractal analyses shows that such approach would lead to the 
same effect as if one may hear, for example, only the first violin of Mahler's music. 

So, it is important to understand when the restricted set of dynamical variables will admit to 
describe the process (state) completely. The same problem was solved in statistical physics, where the 
'rough' description by a restricted number of (thermodynamical) parameters is a basis of its success, 
see the discussion of rough variables description e.g. in the review of Uhlenbeck ||. We will search the 
same solution desiring to build a complete theory of the VHM hadron reactions. 
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We want to note that just VHM process may be in this sense 'simple': at all evidence, the system 
becomes 'quiet' in the VHM region and by this reason its 'rough' thermodynamical description is 
available. It seems natural, therefore, to start investigation of multiple production phenomena from 
the (extremely rare) VHM processes. 



2.1.2 Soft channel of hadron production 

The dominant inelastic hadron processes at n ~ n(s) are saturated by production of low transverse 
momentum hadrons pTj| . One of the approaches explains this phenomena by the nonperturbative effect 
of quarks created from vacuum. 

Corresponding dynamics looks as follows. At the expanse of transverse kinetic motion color charges 
may separated at large distances. Nevertheless the transverse motion is suppressed since separation 
leads to increasing polarization of vacuum, because of confinement phenomenon. Then, as in QED 
JO], the vacuum becomes unstable in regard to the tunneling creation of real fermions. Just on their 
creation, the transverse kinetic energy is spent and, as a result, particles cannot have, with exponential 
accuracy, high transverse energies. This picture is attractive being simple and transparent, but despite 



numerous efforts 12 , there is no quantitative description of this phenomenon till now. Briefly, the 
problem is connected with the unknown mechanism of strong coloured electric fields formation among 
distant colored charges, see also 

One may use other terms. The soft channel of multiple production means the long-range correlation 
among hadron coloured constituents. Under this special correlation the non-Abelian gauge field theory 
conservation laws constraint were implied. They are important in dynamics since each conservation 
law decreases the number of independent degrees of freedom at least on one unity (this may explain 
why hadrons n(s) << n max (s)) 1 i.e. it has the nonperturbative effect. Moreover, in so-called integrable 
systems each independent integral of motion (in involution) reduces the number of degrees of freedom 
on two units. As a result there is no thermalization in such systems ^3 and the corresponding mean 
multiplicity n(s) should be equal to zero, see Appendix K. 

Existence of multiple production, n(s) >> 1, testifies the statement that the thermalization phe- 
nomena exist in hadron processes, i.e. the system of Yang-Mills fields is not completely integrable. But 
the most probable process with n ~ n(s) did not lead to the final state with maximal entropy since 
n(s) << n max , i.e. the definite restrictions on the dissipation dynamics should be taken into account. 
Such problems, being intermediate, are mostly complicated ones. 

The quantitative theory of this phenomena may lead to deep revision of the main notions of existing 
quantum field theory jlM |l6| , see Appendix K. So, the dynamical display of hidden conservation 
laws of the hadron system are probably unstable since we expect that the system is not completely 
integrable, solitary field configurations u c (x,t) fl^| . Then the quantum theory should be able to 
describe quantum excitations of these fields, i.e. to count the fluctuations of 'curved' manifolds. The 
canonical perturbation theory methods, formulated in terms of creation and annihilation of particles 
in the external field u c (x, t), are too complicated, see fis)] and references cited therein. For this reason, 
existing calculations usually do not exceed the semiclassical approximation. We hope that, as described 
in Appendix K, the quantization scheme would be able to solve this problem (see also the example 
described there). 

Another approach assumes that the special 'i-channel' ladder type Feynman diagrams are able to 
describe the n ~ n(s) region |is|j . This approach did not take into account the confinement nonper- 
turbative effects introducing the hardly controllable supposition that the free quarks and gluons may 
form a complete set of states. Formally this is right, but at all evidence, the decomposition on this 
Fock basis is realized in the non-Abelian gauge field theories on zero measure p6l ||cj (see Appendix 
K). Nevertheless one may reject this argument assuming that the process is happening on a sufficiently 
small distances. 

Corresponding contribution came from the so called 'hard Pomeron' ^J. But the intrinsic problems 
of the accuracy of chosen logarithmic approximations p2[ , the understanding of the so called non- 
logarithmic corrections (23[, of the fate of the infrared divergences stay unsolved till now. 
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2.1.3 Multiple production as a process of dissipation 

So, the multiple production of soft hadron phenomena seems unsolvable on the day-to-day level of 
understanding and we may with pure conscience move it away. This is why the hadron inelastic 
reactions lost some popularity, migrating the last two decades to the class of 'non-interesting' problems. 
Yet, in a number of modern fundamental experiments, multiple production plays, at least, the role of 
background to the investigated phenomena and for this reason we should be ready for the quantitative 
estimation of it. 

Our hope to describe such a complicated problem as the multiple production phenomenon in the 
VHM conditions is based on the following idea. At the very beginning of this century, a couple P. 
and T.Euhrenfest, had offered a model to visualize Boltzmann's interpretation of the irreversibility 
phenomena in statistics. The model is extremely simple and fruitful p4fl . It considers the two boxes 
with 2iVb numerated balls. Choosing the label of the balls randomly one must take the ball with the 
corresponding label from one box and put it into another one. One may repeat this action an arbitrary 
number of times t. 

Starting from the highly 'nonequilibrium' state with all balls in one box, Nf, >> 1, it is seen 
stationary with t tendency to equalization of the number of balls in the boxes Fig.l. The stationarity 
means that the number of balls in the other box rises ~ t at least on an early stage of the process. This 
signifies presence of an irreversibly flow (of balls) toward the preferable (equilibrium) state. One can 
hope pl| that this model reflects a physical reality of nonequilibrium processes with initial state very 
far from equilibrium. A theory of such processes with (irreversible) flow toward a state with maximal 
entropy should be sufficiently simple being close to the stationary Markovian. 

The VHM production process may be, at least on an early stage, stationary Markovian. If this is 
so then one may neglect long-range effects, non-perturbative as well, since they are not Markovian as 
follows from the experience described in Sec. 2. 1.2. 

This is possible if the VHM process is happening so fast (being the short-range phenomenon) 
that the confinement forces became 'frozen'. It can be shown that the quantitative reason for this 
phenomena is a fast (exponential) reduction with n of the soft channel contribution into the hadron 
production process. So, we expect a change of the multiple production dynamics in the VHM region. 

Thus, the main input idea consists from two general propositions. The first of them is following: 

(I) The hadron VHM production processes should be close to the stationary Markovian. 
'Freezing' the confinement constraints, the entropy S may exceed for given energy y/s its available 

maximum in the VHM domain. Then one can assume that the VHM final state is in 'equilibrium', or 
is close to it. So, 

(II) The VHM final state should be close to equilibrium 
is our second basic proposition. 

We would select and appreciate particle physics models in accordance with these propositions. 



2.2 S- matrix interpretation of thermodynamics 

The field-theoretical description of statistical systems at a finite temperature is based usually on the 
formal analogy between imaginary time and inverse temperature (3 = 1 /T p6[ . This analogy is formu- 
lated by Schwinger (l|] as the 'euclidean postulate' and it assumes that (i) the system is in equilibrium, 
i.e. it should allow the arbitrary rearrangement of states of temporal sequence in the described pro- 
cess^], and (ii) there is not special space-time long-range correlations among states of the process, i.e, 
for instance, the symmetry constraints should not play a crucial role. We do not know ad hoc whether 
or not to apply the 'euclidean postulate' for given n and s, even if ( |l.l| ) is satisfied. For this reason we 
are forced to formulate the theory in the natural real-time terms. 

The first important quantitative attempt to build the real-time finite-temperature field theory [ p7| 
discovered the formal problem of the so called 'pinch-singularities'. Further investigation of the theory 

2 One can say that the opposite flow is never seen. 'What never? No never! What never? Well, 
hardly ever.' This dialog was taken from |2[| 

3 In other terms, one may have the possibility to apply the ergodic hypothesis. 
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has allowed to demonstrate the cancellation mechanism of these unphysical singularities p8| . This 
is attained by doubling the degrees of freedom |^9|, 30 : the Green functions of the theory represent 
2x2 matrix plj . It surely makes the theory more complicated, but the operator formalism of the 
thermo-field dynamics |32[] shows the unavoidable character of this complication. 



Notice that the canonical real-time finite temperature field-theoretical description |29, M of the 



statistical systems based on the Kubo-Martin-Schwinger (KMS) [|29|, |33j, |34j boundary condition for a 
field is: 

$(*) = $(* - i0). (2.2) 

It, without fail, leads to the equilibrium fluctuation-dissipation conditions plj ] (see also [p6|). By this 
reason it can not be applied in our case, where the dissipation problem is solved. The origin of this 
boundary condition is shown in Appendix A. 

We will use a more natural microcanonical formalism for particle physics^. The thermodynam- 
ical 'rough' variables are introduced in these approach as the Lagrange multipliers of corresponding 
conservation laws. The physical meaning of these 'rough' variables are defined by the corresponding 
equations of state. 

We shall use the iS-matrix approach which is natural for the description of the time evolution. (The 
S'-matrix description was used also in (3^, pjj.) For this purpose the amplitudes 

^ Pi i P2 , • ■ ■ , Pn \ qii 92 , ■••) Qm = Q"nm (P,q) (2-3) 

of the to- into n-particles transition will be introduced. The in- and out-states must be composed from 
mass-shell particles [|o| . Using these amplitudes we will calculate 

R nm (p,q) = \a nm (p, q)\ 2 = 
=< Pi,P2, ...,p n \qi,q2, -,q m >< 91,92, ...,q m \pi,P2, -,Pn > ■ (2.4) 

This will lead to the doubling of the degrees of freedom. 

The temperature description will be introduced (see also ]4l| ) noting that, for instance, 

dT n = \a nm {p,q)\ 2 d£l n (p), 

dVnip) = n" vAty e ^ = ( p2 + m2 ) 1/2 ' ( 2 - 5 ) 

is the differential measure of the final state. It is a first example where the usefulness of the probability- 
like quantity ~ |a nm | 2 is seen. 

The measure (B7q) is defined on the energy-momentum shell: 



5>i = P. (2.6) 

i=l 

It should be underlined that a nm (p, q) are the translationally invariant amplitudes and four equalities: 

n m 



are obeyed identically. So, (2.6) are the constraints and to take them into account one may multiply 
dT n on 

n 

fe=i 

where a is the time-like 4-vector. It is evident that integration over a with factor e~ %aF gives the 
constraints (|2.6|). 



The statistical methods for particles physics was discussed also in WA 
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One may simplify the calculation assuming that all calculations are performed, for example, in the 
CM frame P = (E,0). Then one may ignore the space components considering a — (ao,0). This is 
the equivalent of the assumption that only the energy conservation law is important. 

The last step is the substitution a® = if3, where /3 is our Lagrange multiplier. To define its physical 
meaning one should solve the equation of state: 

E = -j8 ln I dTn n e_ " £(pfc) = -^ in / dr ^- ( 2 - 8 ) 

k—l 

Such a definition of temperature as the Lagrange multiplier of the energy conservation law is obvious 
for microcanonical description fl33| . 

The initial-state temperature will be introduced in the same way, taking into account fl2.7| ). So, 
we will construct the two-temperature theory. It is impossible to use the KMS boundary condition 
in such a two temperatures description (the equation of state can be applied at the very end of the 
calculations) . 

It should be notice that the 'density matrix' R n , m (p, <?), defined in ( [2.4|) , describes the 'closed-path 
motion' in the functional space. So, if 

< pi,p 2 , —,p n \qi,Q2, ■•■,9m >=< n, out\e' s ^ +) \m, in > (2.9) 

and 

< Pl,P2, ...,Pn\qi,q2, ■■■,9m >* = < <7l,<?2, ., q m \Pl , P2, -,Pn > = 

= <m,in\e~ lS{<s '- ) \n,out> (2.10) 

then, by definition, 

$+(a DO ) = $_(a DO ) = $( C r 0O ), (2.11) 

with some 'turning-point' fields < &((t 00 ), where is the remote hypersurface. The value of ^((Too) 



specifies the environment of the system. We will show that (2.11) coincides with the KMS boundary 



condition in some special cases. Here consequences of the vacuum boundary condition: 

<S>(cr x )= (2.12) 

are analyzed. 

One should admit also that below boundary conditions are not unique: one can consider arbitrary 
organization of the environment of the considered system. The 5-matrix interpretation is able to show 
the way as an arbitrary boundary condition may be adopted. This should extend the potentialities of 
the real-time finite-temperature field-theoretical methods. 

2.2.1 Example 

It seems useful to illustrate the above microcanonical approach by the simplest example, see also @. 
By definition, the n particles production cross section 

/n 
dQ n {p)S{qi +q 2 - ^2pi)\a n (p, q)\ 2 , (2.13) 

i—l 

where a n (p, q) = &„ 2 (Pj ?) is the ordinary n particle production amplitude in accelerator experiments. 

Considering the Fourier transform of energy-momentum conservation <5-function one can introduce 
the generating function p n , see |41 and references cited therein^. We may find in the result that a n is 



defined by the equality 



/-|-'iOO j o 
E = e( qi )+e(q 2 ), (2.14) 
-too 27r 



3 The generating functionals method was developed in [|42| 
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where 



Pn(P) 




I 2 = / dT n ((3). 



(2.15) 



Most probable value of (3 in (2.14) is defined by the equation of state (2.8). Inserting ( 2.1 5| ) into ( |2.14 ) 
we find the expression (2.13) if the momentum conservation shell is neglected. The last one is possible 



since the cross sections are always measured in the definite frame. 
Let us consider the simplest example of noninteracting particles 

Pn (p) = {2nmK 1 (Pm)/p} n , 



11: 



where K\ is the Bessel function. Inserting this expression into (2.8) we can find that in the nonrela- 
tivistic case (n ~ n max ) 

3 (n-1) 



Pc 



i.e., we find the well known equality: 



2 (y/s — nm) ' 



E kin - 2 T, 



(2.16) 



where Ekin = {\fs—nm)/{n—l) is the mean kinetic energy and T = l/p c is the temperature (Boltzmann 
constant was taken equal to one). 



It is important to note that the equation (2.8) has a unique real solution p c (s, n) rising with n and 
decreasing with s |33] . 



The expansion of the integral (2.14) near /3 c (s, n) unavoidably gives an asymptotic series with zero 
convergence radii since p n (P) is the essentially nonlinear function of /?, see also Sec. 2. 2. 2. This means 
that, generally speaking, fluctuations in the vicinity of f3 c (s,n) may be arbitrarily high and in this 
case (3 c (s,n) has no physical sense. But if the fluctuations are Gaussian, then p n (P) coincides with the 
partition function of the n particle state and j3 c (s,n) may be interpreted as the inverse temperature. 
We will put the observation of this important fact in the basis of our thermodynamical description of 
the VHM region. 



2.2.2 Relaxation of correlations 

The notion of 'equilibrium' over some parameter X in our understanding is a requirement that the 
fluctuations in the vicinity of its mean value, X, have a Gaussian character. Notice, in this case, one 
can use this variable for a 'rough' description of the system. We would like to show now that the 
corresponding equilibrium condition would have the meaning of the correlations relaxation condition 
of Bogolyubov |42]fl, see also . Let us define the con dition s when the fluctuations in the vicinity of 
Pc are Gaussian |44|. Firstly, to estimate the integral ( 2.14 ) in the vicinity of the extremum, p c , we 
should expand In p n (/3 + f3 c ) over (3: 



1 d 2 



if 



and, secondly, expand the exponent in the integral ( [2.14 ) over, for instance, 

a 3 ln Pll (/3 c )/c)/3 c 3 ,..., 



(2.17) 



etc. In the result, if higher terms in ( p.l7| ) are neglected, k-th term of the perturbation series 

k 



Pn,k 



f d 3 \np n (p c )/dp 3 c 
\(dHnp n (p c )/dp2)a/z 



3k + 1 



(2.18) 



3 The term 'vanishing of correlations' was used by N.N. Bogolyubov for this phenomena. 
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Therefore, because of Euler's T((3fc + l)/2) function, the perturbation theory near j3 c leads to the 
asymptotic series. The supposition to define this series formally, for instance, in the Borcl sense is not 
interesting from physical point of view. Indeed, such a formal solution assumes that the fluctuations 
near C may be arbitrarily high. Then, by this reason, the value of (3 C loses its significance: arbitrary 
values of ((3 — j3 c ) are important in this case. 

Nevertheless it is important to know that our asymptotic series exists in some definite sense, i.e. 
we can calculate the integral over (3 by expanding it over (f3 — (3 C ). Therefore, if the considered series 
is asymptotic, we may estimate it by first term if 

d 3 \n Pn (l3 c )/dl3 3 c « (a 2 lnp„(/3 c )/9/3 c 2 ) 3 / 2 . (2.19) 

One of the possible solutions of this condition is 

d 3 ln Pn {f3 c )/dl3 3 c ^0. (2.20) 

If this condition is satisfied, then the fluctuations are Gaussian with dispersion 

-{c^lnp^)/^ 2 } 1 / 2 , 



see ( p.lTp . 

Let us consider now (2.2C) carefully. We will find by computing derivatives that this condition 
means the following approximate equality: 

_L_ 3 _^_ + 2 M« , (2.21) 
Pn Pi Pt 

(k) 

where p n means the k-th derivative. For identical particles, 

pW(J3 c ) = n k (-l) k J dT n (/3 c ) IJefa), (2-22) 



Therefore, the left hand side of ( 2.21 ) is the 3-point correlator K% since <fT„(/3 c ) is a density of states 
for given j3: 



K 3 = Jdn 3 (q) (< UU<Qi) >A, -3 < Iltiefe) >/3c< c(9s) >/3 c + 

+ 2llL<<ft)>/»o). (2-23) 



where the index (3 C means that averaging is performed with the Boltzmann factor exp{— (3 c e(q)}. 

Notice, in distinction with Bogolyubov, Kg is the energy correlation function. So, in our interpre- 
tation, one can introduce the notion of temperature l//? c if and only if the macroscopic energy flows, 
measured by the corresponding correlation functions, are to die out. 

As a result, to have all the fluctuations in the vicinity of j3 c Gaussian, we should have K m « 0, 
m > 3. Notice, as follows from ( 2.19| ), the set of minimal conditions actually looks as follows: 

K m « K 2 , m > 3. (2.24) 

If the experiment confirms this conditions then, independently from the number of produced particles, 
the final state may be described with high enough accuracy by one parameter (3 C and the energy 
spectrum of particles is Gaussian. In this conditions one may return to the statistical |Q and the 
hydrodynamical models [45fl . 

Considering f3 c as a physical (measurable) quantity, we are forced to assume that both the total 
energy of the system, y/s = E, and the conjugate to it, variable C , may be measured simultaneously 
with high accuracy. 
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2.2.3 Connection with Matsubara theory 

We would like to show now that the ordinary big partition function of the statistical system coincide 
with 



E 



(2.25) 



(/3i,zi;/32,Z2) 



where R nm (p, q) is defined by ( |2.4| ) . The summation and integration are performed with constraints 
that the mean energy of particles in the initial(final) state is l//3i(l//?2). On may interpret 1//3 in the 
first approximation as the temperature and Z\{zn) as the activity for initial(final) state. 

Direct calculation, see Appendix B, gives the following expression for generating functional: 



p(p,z) = e-™^Ro(<l>), 
where the particle number operator (4>(x) = 8/5<f>(x)) 



and 



- J dxdx' '(</> + (x)-D_| (x — x', (3%, ^)0- {x') 

~4>-(x)D- + (x - x',(3i, zi)4> + {x')) 

R (<!>)=Z(<j> + )Z*(-<f>^, 



(2.26) 

(2.27) 
(2.28) 



where Z(cj>) is defined in ( B.10 ): 



z{<j>) = J D ^sm-iv{^) 

and, for the vacuum boundary condition $((7oo) = 0, 

D + _(x-x',l3,z) = -i J dD, 1 {q)e iq{ - x - x 'h- <' l h{q) 
D- + (x-x',/3,z)=i [ dnx{q)e- iq ^- x '^e-^z{q) 



(2.29) 
(2.30) 



are respectively the positive and negative frequency correlation functions at z = 1. 
It is evident, 



" f 

R nm (p,q) = J] e ft 

k=l 



Szi(j>k) 



n 

fc=i 



5z 2 {qk) 



p(fl,z) 



(2.31) 



Zi=0 



Notice, defining R n m{p,q) through the generating functional we extract the Boltzmann factors e _/3e 
since the energy-momentum conservation 5-functions were extracted from amplitudes a nm (p,q). 

We suppose that Z{<t>) may be computed perturbatively. As a result, (J — 5/Sj is the variational 
derivative) 

R(P Z ) — e - iVr (-y+)+ iy (-y-) e 4 / dxdx 'ji(x)D ik (x-x':l3,z)j k (x') ^ 

where Dn~(x — x') is the matrix Green function. These Green functions are defined on the Mills p6| 
time contours C± in the complex time plane (C_ = C*), see Fig. 2. This definition of the time contours 
coincides with the Keldysh' tim e con tour [p0| . 

The generating functional (|2.32[ ) has the same structure as the generating functional of Niemi 
and Semenoff [P8| . The difference is on ly in the de finitio n of Green functions D^. This choice is a 
conseq uence of the boundary condition ( |B.6| ). So, if ( B.32 ) is used, then the Green function is defined 
by eq.( B.45 ). Notice also that if j3\ = 02 = /3 then a new Green function obey KMS boundary condition, 
see ( |B.4SD . 

Following Niemi and Semenoff |2gj one can write ( 2.32| ) in the form: 



PC8)= / D NS ^ S «^\ 



(2.33) 
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where the functional measure Dns<& and the action Sns(^) are defined on the closed complex time 
contour CnSi see Fig. 3. The choice of initial time U and tf is arbitrary. Then one can perform shifts: 
ti — ► — oo and tf —* +00. In result, (i) if Pi — P2 — 0, (ii) if contributions from imaginary parts 



C-\ and C-p of the contour C'ns have disappeared in this limit, (iii) if the integral (2.33) may be 

calculated perturbatively then this integral is a compact form of the representation ( 2.32] ). 



Notice that the requirements (i) - (iii) are the equivalent of the Euclidean postulate of Schwinger. 
In this frame one can consider another limit tf — > ij. Then the Cns contour reduces to the Matsubara 
imaginary time contour, Fig. 4. 

Later on we will use this S'-matrix interpretation of thermodynamics. But one should take in mind 
that corresponding results will hide assumptions (i) - (iii) . 

We would like to mention the ambivalent role of external particles in our S-matrix interpretation 
of thermodynamics. In the ordinary Matsubara formalism the temperature is measured assuming 
that the system under consideration is in equilibrium with the thermostat, i.e. temperature is the 
energy characteristics of interacting particles. In our definition the temperature is the mean energy of 
produced, i.e. non-interacting, particles. It can be shown that both definitions lead to the same result. 

Explanation of this coincidence is the following. Let us consider the point of particle production 
as the coordinate of fictitious 'particle'. This 'particle' interacts since the connected contributions into 
the amplitudes a nm only are considered, and has the equal to produced particles momentum and so 
on. The set of these 'particles' form a system. Interaction among these 'particles' may be described 
by corresponding correlation functions, see Sec. 2. 3. 3. 



Let us consider now the limit tf — > ti. In this limit (2.33) reduces to 

p(f3) = J D M $e- s **W, (2.34) 
where the imaginary time measure Dm® and action Sm(®) are defined on the M atsu bara time contour. 



The periodic boundary condition (2^2) shou ld be used calculating integral ( 2.34 ). The rules and 



corresponding problems in the integral ( 2.34 ) can be calculated as described in many textbooks, see 
also WA. 



In the limit considered the time was eliminated in the formalism and the integral in ( 2.34 ) performed 
over all states of the 'particles' system with the weight e^ Su ^\ Notice, the doubling of degrees of 
freedom has disappeared and our fictitious 'particles' became real ones. 

On the other hand, the produced particles may be considered as the probes through which we 
measure the interacting fields. As was mentioned above, their mean energy defines the temperature, 
if the energy correlations are relaxed. If even one of the conditions (i) - (iii) is not satisfied then one 
cannot reduce our S'-matrix formalism to the imaginary time Matsubara theory. Then one can ask: is 
there any possibility, staying in the frame of S'-matrix formalism, to conserve the statistics formalism. 
This question is discussed in the Appendix C, where the Wigner functions approach is applied. It may 
be shown that the formalism may be generalized to describe the kinetic phase of the nonequilibrium 
process, where the temperature should have the local meaning E^j. The comparison with the 'local 
equilibrium hypothesis' is discussed at the end of Appendix C. 

2.3 Classification of asymptotics over multiplicity 

Our further consideration will based on the model independent (formal) classification of asymptotics 
& 

2.3.1 'Thermodynamical' limit 

We will consider the generating function: 

Umax 

T(s, z) = ^2 z™er„(s), s=(pi+ p 2 ) 2 » to 2 , n max = y/s/m. (2.35) 
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This step is natural since the number of particles is not conserved in our problem. So, the total cross 
section and the averaged multiplicity will be: 



vtot(s) = T(s, 1) = er n (s), n(s) 



n(a n (s) / a to t) 



dz 



lnT(s, z) 



z = l 



At the same time, the inverse Mellin transform gives 



1 d n 

n\ dz n 



T{s,z) 



z=0 



\ ( dz , , 1 
77- <P — rr T (s,z) = — 



2iri 



^£ e (— n In 2+ln T(s,z)) 



The essential values of z in this integral are defined by the equation (of state): 

n = z— hxTiz, s). 
oz 



(2.36) 



(2.37) 



(2.38) 



Taking into account the definition of the mean multiplicity n(s), given in (2.36), we can conclude that 
the solution of (2.38) z c is equal to one at n = n(s). Therefore, z > 1 are essential in the VHM domain. 

The asymptotics over n (n << n max is assumed) are g overned by the smallest solution z c of (2.38 ) 
because of the asymptotic estimation of the integral ( [2.37 ) : 



a n (s) oc e -«l»*=Ks). 



(2.39) 



Let us assume that in the VHM region and at high energies, y/s — > oo, there exist such a value of 
z c (n, s) that we can neglect in ( [2.35 ) the dependence on the upper boundary n max . This formal trick 
with the thermodynamical limi t allows to consider T(z, s) as the nontrivial function of z for finite s. 

Then, it follows from (|2.38|) that 



z c (n, s) — » z s at n e VHM, 



(2.40) 



where z s is the leftmost singularity of T(z, s) in the right half plane of complex z. One can say that 
the singularity of T(z, s) attracts z c (n, s) if n £VHM. We will put this observation in the basis of VHM 
processes phenomenology. 

We would like to underline once more that actually T(z,s) is regular for arbitrary finite z if s is 
finite. But z c (n, s) behaves in the VHM domain as if it is attracted by the (imaginary) singularity z s . 
And just this z c (n, s) defines a n in the VHM domain. We want to note that actually the energy s/s 
should be high enough to use such an estimation. 



2.3.2 Classes and their physical content 



One can notice from the estimation ( |2.39 ) that er„ weakly depends on the character of the singularity. 



Therefore it is enough to classify only the possible positions of z s . We may distinguish following 
possibilities: 

(A) z s = oo : er„ < 0(e~ n ) 

(B) z s = l: a n > 0(e- n ) 

(C) 1< z s < oo : a n = 0{e~ n ), (2.41) 

i.e., following this classification, the cross section may decrease faster (A), slower (B), or as (C) an 
arbitrary power of e~ n . It is evident, if all these possibilities may be realized in nature, then we should 
expect the asymptotics (B). 

As was explained in Sec. 2. 2.1, a n has the meaning of the n particle partition function in the energy 
representation. Then T(z, s) should be the 'big partition function'. Taking this interpretation into 
account, as follows from Lee- Yang theorem |5l|], T(z,s) can not be singular at \z\ < 1. 

At the same time, the direct calculations based on the physically acceptable interaction potentials 
give the following restriction from above: 

(D) a n < 0(l/n) (2.42) 
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This means that a n should decrease faster than any power of 1/n. 

It should be noted that our classification predicts rough (asymptotic) behavior only and did not 
exclude local increase of the cross section <j n . 

One may notice that 

_ l h Mi =lnz c (n,s) + 0(l/n ). (2.43) 
n <T t ot{s) 

Using thermodynamical terminology, the asymptotics of a n is governed by the physical value of the 
activity z c (n, s). One can introduce also the chemical potential /i c (n, s). It defines the work needed for 
one particle creation, lnz c (n, s) = f3 c (n, s)ii c {n, s), where e(n, s) = 1/ (3 c (n, s) is the produced particles 
mean energy. So, one may introduce the chemical potential if and only if f3 c (n, s) and z c (n, s) may be 
used as the 'rough' variables. 

Then the above formulated classification has a natural explanation. So, (A) means that the system 
is stable with reference to particle production and the activity z c {n, s) is the increasing function of n, 
the asymptotics (B) may realized if and only if the system is unstable. In this case z c (n,s) is the 
decreasing function. The asymptotics (C) is not realized in equilibrium thermodynamics p2[ . 

We will show that the asymptotics (A) reflects the multiperipheral processes kinematics: created 
particles form jets moving in the CM frame with different velocities along the incoming particles di- 
rections, i.e. with restricted transverse momentum, see Sec. 3. 1.1. The asymptotics (B) assumes the 
condensation-like phenomena, see Sec. 3. 3. The third type asymptotics (C) is predicted by station- 
ary Markovian processes with the pQCD jets kinematics, see Sec. 3. 2. 2. The DIS kinematics may be 
considered as the intermediate, see Sec. 3. 2.1. 



This interpretation of classes ( 2.41 ) allows to conclude that we should expect reorganization of 
production dynamics in the VHM region: the soft channel (A) of particle production should yield a 
place to the hard dynamics (C), if the ground state of the investigated system is stable with reference 
to the particle production. Otherwise we will have asymptotics (B). 

2.3.3 Group decomposition 

Let us consider the system with several correlation scales. For example, in statistics one should dis- 
tinguish correlation length among particles (molecules) and correlation length among droplets if the 
two-phase region is considered. In particle physics, one should distinguish in this sense correlation 
among particles produced in result of resonance decay and correlations among resonances. In pQCD 
one may distinguish correlations of particles in jet and correlation among jets. 

There exist many model description of this physical picture. In statistics Mayer's group decom- 
position |^5| is well known. In particle physics one should note also the many-component formalism 
We will consider the generating functions (functionals) formalism g| considering mostly jet 
correlations. In many respects it overlaps the above mentioned approaches. 

The generating function T(z, s) may be written in the form: 

00 / 1 \k 00 

lnT(z, s)=J2 ^-^p-Cfe( s ) = E ^ ( 2 - 44 ) 

k=l ' l=l 

where the coefficients C k are the moments of the multiplicity distribution 

P n {s) = <Tn{s)/a to t{s). (2.45) 

So, 

d{s) =J2nP n (s) =n(s), C 2 (s) = ^2n(n-l)P n (s)-n(sf (2.46) 

n n 

and so on. Using the connection with the inclusive distribution functions fk(Qi,Q2, —,Qk)'- 

00 (z - i) k r 

T(z,s) = J2 fcj dn k (q)f k (q u q 2 ,...,q k ;s), (2.47) 

k=l ' ^ 

7 The example considered in Etq] illustrate this approach. 
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it is easy to find that 

Ci(*)=/tff2 1 (?)/i(?;a)=/i(«), 
C 2 (s) =Jdn 2 (q){f 2 ( qi ,q 2 ;s) - Mqv, «)/i(<ft; a)} = h(s) - f?(s), (2.48) 

etc. Generally, 



^*(-)=E^ £ ME*<-* n{^}- ( 2 - 49 ) 

(=1 {fe}i=0 \i=l / i=l k J 

where {fc}/ = k\, k 2 , fcj and 

fk(s) = J dQ k (q)f k (q 1 ,q 2 , ...,q k ;s). 

One may invert formulae ( |2.4S| ): 

{n fc }=0 \fc=l / k=l V 7 



The Mayer's group coefficients bi in (2.44) have the following connection with Ct- 



k=0 

It seems useful to illustrate the effectiveness of the generating function method by the following 
example. We will consider the transformation (multiplicity n — > activity z) to show the origin of the 
Koba-Nielsen-Olesen scaling (KNO-scaling)^. 

If C m — 0, m > 1, then <j n is described by the Poisson formulae: 

a n (s)=a tot (s)e- n1 ^-. (2.52) 

It corresponds to the case of absence of correlations. 
Let us consider more weak assumption: 

C m (s) = lm (Ci(s)) m , (2.53) 

where 7 m is the energy independent constant, see also [53] ], where a generalization of KNO scaling on 
the semi-inclusive processes was offered. Then 



lnT(z,s)=^2^{(z-l)n(s)r, (2.54) 



To find the consequences of this assumption, let us find the most probable values of z. The equation 
of state 

n = z— hiT(z, s) 
oz 

has solution z(n, s) increasing with n since T(z, s) is an increasing function of z, if and only if, T(z, s) 
is nonsingular at finite z. As was mentioned above, the last condition has deep physical meaning and 
practically assumes the absence of the first order phase transition j51[ . 
Let us introduce a new variable: 

A = (z - l)n(s). (2.55) 



8 In private discussion with one of the authors (A.S.) in the summer of 1973, Z.Koba noted that the 
main reason of investigation leading to the KNO-scaling was just the generating functional method of 
Bogolyubov [E§ 
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The corresponding equation of state looks as follows: 

^ T =fl+^ T N )^-lnT'(A). (2.56) 
n(s) \ n{s) J o\ 

So, with 0(X/n(s)) accuracy, one can assume that 

A ~ \ c (n/h(s)). (2.57) 

are essential. It follows from this estimation that such scaling dependence is rightful at least in the 
neighborhood of z = 1, i.e. in vicinity of main contributions into Otot- This gives: 

n(s)a n (s) = a to t{s)^{n/n(s)), (2.58) 

where 

ip(n/n(s)) ~ T{\ c (n/n(s))) exp{n/n(s)A c (n/n(s))} < 0(e~ n ) (2.59) 

is the unknown function. The asymptotic estimation follows from the fact that A c = \ c {n/n{s)) should 
be a nondecreasing function of n, as follows from nonsingularity of T{z, s). 

The estimation ( 2.57 ) is right at least at s — > oo. The range validity of n, where solution of ( 2.57 ) 
is acceptable, depends from exact form of T(z,s). Indeed, if lnT(z) ~ exp{7A(z)}, 7 = const > 0, 
then (2.57) is right at all values of n and it is enough to have the condition s — > 00. But if lnT(z, s) ~ 
(1 + aA(z)) 7 , 7 = const > 0, then (2.57) is acceptable if and only if n « n 2 (s). 

Representation (2.58) shows that just n(s) is the natural scale of multiplicity n [g4J. This repre- 
sentation was offered first as a reaction on the so called Feynman scaling for inclusive cross section: 

k 

Ate, <&,..., ft) ~n -rv ( 2 - e °) 

.=1 e[ - qi > 



As follows from estimation (2.59), the limiting KNO prediction assumes that a n = 0(e n ). In this 
regime T(z, s) should be singular at z = z c {s) > 1. The normalization condition 

dT(z, s) _ 
—^\ z ^ = n(s) 

gives: z c {s) = 1 + 7/n(s), where 7 > is the constant. Notice, such behavior of the big partition 
function T(z, s) is natural for stationary Markovian processes described by logistic equations [ ^5| . In 
the field theory such equation describe s the QCD jets [f56| . 

We wish to generalize expansion (2.44) to take into account the possibility of many-component 
structure of the multiple production processes Q. Let us consider particle production through the 
generation, for instance, of jets. In this case decay of a particle of high virtuality \q\ » m forms a 
jet of lower virtuality particles. It is evident that one should distinguish correlation among particles in 
the jet, and correlation among jets. 

Let uj ni {m,i) be the probability that i-ih jet of mass nrii includes rij particles, 1 < n,; < n, where 

^n t = n (2.61) 

i=l 

The jets are the result of particles decay. Then let us assume that N 1 (mi,Pi) defines the mean number 
of jets of mass m, and momentum pf. 

X Sllflx {^^e-^N^p^irm)} , (2.62) 

where {n}^. — (711,712, ...,njv, )- Notice the Boltzmann fa ctor e~^ e , where e(p) = m + p 2 /2m is the jets 
energy, play the same role as the corresponding factor in ( [j.l5| ) and introduced to take into account the 
energy conservation law. We consider the VHM domain and for this reason (p 2 /2m) << 1 is assumed. 
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It is useful to avoid the particles number conservation law ( [2.61 ) . For this purpose we will introduce 
(0 W(/3,z)=^p( 1 )(/3)=exp|y' ^^e-^N x {m,p){t{z,m)-l)}, (2.63) 

where 



t(z, to) = z n uj n (m). 



(2.64) 



Comparing (2.63) with ( |2.44 ) we may conclude that t(z,m) plays the role of activity of jets. Then the 
generalization is evident: 

= exp {E fc / Ilti {dmmi{Pi)e-^\t{z, m 4 ) - 1)} x 

xN k (m 1 ,p 1 ,...,m k ,p k )} , (2.65) 

where N k has the same meaning as C k , i.e. N k is the correlation function of k jets. 

2.3.4 Energy- multiplicity asymptotics equivalence 

Let us consider the following 'bootstrap' regime when p(/3, z) is defined by the equation: 



pGM« / ^^e-^N x {m,p)t{z,m). 
J 2m (2iry 



(2.66) 



Inserting here the st rict e xpression (2.65) we find a nonlinear equation for t(z, m) 
The solution of ( 2 . 66| ) assumes that 



dm d 3 p 
2m (2tt) 3 



» 



dm d 3 p , 6(p) 



Nit 



and 



dm d 3 p fl , - 
2m (2?r) 3 



2to (2tt) 3 
dm i d 3 pi -pe t \ ( dm 2 d 3 p 2 n _ p 



2toi (2tt) 



2to 2 (2tt) 3 



e~ p H N- 



(2.67) 



(2.68) 



To solve eq.(2.66) in the VHM region, where the leftmost singularity over z is important, let us 
consider the anzats: 

(1 — [z — l)a(m)) K ° 

where <p(z, m) is the polynomial function of z, ip(z — 1, m) = 1. Using the normalization condition: 

d 



n 3 = ^( 2 > m ) 



d 

a(m)K = rij - if' {I, to), ip'(l,m)= — ip(z,m) 



we can find: 

a(m)n = hj r ^, ..„„ y. ^, — , 

The partition function of the jet m) defined by anzats ( 2.69| ) is singular at 

1 



z=l 



z s {m) = 1 + 



a(m) 



(2.70) 
(2.71) 

(2.72) 



This singularity would be significant in the VHM region if z s (m) is decreasing function of m. This 
means an assumption that 

ip'(l,m) 



nj(m) 



at to — ► oo . 
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So, in first approximation we will choose 



a(m) = rij(m). 

This choice may be confirmed by concrete model calculations. 



(2.73) 



Taking into account the energy conservation law, conditions (2.67) and (2.68) are satisfied if 



exp 



n n j (s)-n j (s/4) \ 
%(s)%(s/4) J 



<< 1. 



(2.74) 



at n GVHM. Therefore, ( 2.69 ) obey eq.( 2.66| ) with exponential accuracy in the VHM region, i.e. if 

fij(s)fij(s/4) k q 



n » 



rij {s)-n j (s/A) z s (s/4) - z s (s) ' 



(2.75) 



We assume here that one can find so large n and s that with expo nent ial accuracy the factors ~ Nk 
did not play an important role. But at low energies the condition ( |l.2| ) is important and the factors 
~ Nk should be taken into account. 

Notice now important consequence of our 'bootstrap' solution: it means that we can leave produc- 
tion of the heavy jets only, if n £ VHM . On the other hand, let us choose n = zofij(s), where zq > 1 
is the function of s, and consider s — > oo. Then the condition (2.75) defines zq: if 



z » 



n s (a/4) 



nj(s)-nj(s/A) 



then we are able to obey the inequalities (2.67) and (2.68). 
The jet mean multiplicity, see Sec. 3. 2. 2, 



Then 



r»j(»/4) 



nj(s) -n 3 (s/A) 



In rij(s) ~ vlns. 



'Ins << rij(s). 



(2.76) 

(2.77) 
(2.78) 



at s — > oo. Therefore, (2.76) may be satisfied outside the VHM domain. 

Let us compare now the solutions of the equation of state. Inserting (2.6E) into ( 2.38[ ) we can find 
for a jet of mass y/s that 

1 1 «o 1 



n nj(s) 
The two-jet contribution of the masses ~ y/s/2 gives: 



K 



= 1 



1 



nj(*/4) 



n 



n 



(2.79) 



(2.80) 



At arbitrary finite energies [z\ — z\) > and, as follows from ( [2.7§| ), they decrease ~ (l/nj(s)Vhis) 
with energy. 

Noting the normalization condition, T(z = 1, s) = <J to t(s), and assuming that the vacuum is stable, 
i.e. er„ < 0(e~ n ), we can conclude that 

- if n G VHM then z s attracts z c , i.e. z s — * z c , and if z c — 1 << 1 then this contributions should 
be significant in cr tot ; 

oo, then z c — 1 << 1, and if n satisfy the inequality (2.75), or if zq satisfy the inequality 



if . 



( 2. 76] ), then considered contributions are significant in crtot- 

It is the (energy-multiplicity) asymptotics equivalence principle. One of the simplest consequences of 
this principle is the prediction that the mean transverse momentum of created particles should increase 
with multiplicity at sufficiently high energies. 

This principle is the consequence of independence of contributions in the VHM domain on the type 
of singularity in the complex z plane and of the energy conservation law. Just the last one shifts the 
two-jet singularity to the right and z c (s) < z c {s/A). 
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We would like to notice also that this effect, when the mostly 'energetic population' survives was 
described mathematically by V.Volterra |i5| . Intuitively evident is that one may find the 'energetic 
population' searching the VHM one, or, it is the same, giving it a rich supply, i.e. to give the population 
enough energy. This is our (energy- multiplicity) equivalence ((e — n)-equivalence) principle. 

Notice, if the amount of supply is too high then few populations may grow. This is the case when 
the difference {z 2 c — z\) > tends to zero at high energies. 

We would like to note that singular at finite z partition functions was predicted in the (A0 3 )g- 



theory p7y, in QCD jets pa], in the generalized Bose-Einstein distribution model 58 . In all of this 



models decay of the essentially nonequilibrium initial state (highly virtual parton, heavy resonance, 
etc) was described. 

One may distinguish the phases of the media by a characteristic correlation length. Then the phase 
transition may be considered as the process of changing correlation length. Our 'bootstrap' solution 
predicts just such phenomena: at low multiplicities the long-range correlations among light jets is 
dominant. The 'bootstrap' solution predicts that for VHM processes just the short-range correlations 
among particles of the heavy jet become dominant. The (e — n)-equivalence means that this transition 
is a pure dynamical effect. 



3 Model predictions 

Multiple production phenomena was first observed more than seventy years ago p39[| . During this time a 
vast experimental information was accumulated concerning hadron inelastic interactions, see the review 
papers ]Tc| ]. 

Now we know that at high energy \fs: 

(i) The total cross section o~tot(s) of hadron interactions is enhanced almost completely by the inelastic 
channels; 

(ii) The mean multiplicity of produced hadrons n(s) slowly (logarithmically) grows with i/s; 

(iii) The interaction radii of hadrons b slowly (logarithmically) increase with energy; 

(iv) The multiplicity distribution a n {s) is wider than the Poisson distribution.; 

(v) The mean value of transverse momentum k of produced hadrons is restricted and is independent 
of the incident energy y/s and produced particle multiplicity n; 

(vi) The one-particle energy spectrum da ~ de/e. 

First of all, (i) means that the high energy hadron interaction may be considered as the ordinary 
dissipation process. In this process the kinetic energy of incident particles is spent in produced particle 
mass formation. 

The VHM process takes place in the vacuum and then it was assumed on the early stages that 
the multiple production phenomena reflects a natural tendency of the excited hadron system to get 
to equilibrium with the environment jjj. In this way one can introduce as a first approximation the 
model that the excited hadron system evolves without any restrictions. In this model we should have 
n(s) ~ y/s. The dissipation is maximal in this case and the entropy S exceeds its maximum. This 
simple model has definite popularity up to 70-th. But the experimental data (ii) and (iii) prohibit this 
model and it was forgotten. 

Choosing the model we would like to hope that the considered model 
-takes into account experimental conditions (i) - (vi) in the n ~ n(s) domain; 

- has natural asymptotics over multiplicity to the VHM region. 
It is necessary to remember also that 

- New channels of hadron production may arise in the VHM region. 

It is impossible to understand all possibilities without those offered in Sec. 2. 3. 2 in the classification of 
asymptotics. 

Thus, we will observe predictions of 

- Multiperipheral models, 
distinguishing the soft Pomeron models, see p(| . 

- The dual-resonance model 

predictions for the VHM region are described also. It can be shown that this models predict asymptotics 
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(A) if n >> n(s) 2 . Just this result explains why VHM domain is defined by the condition 

We are forced since it allows to include pQCD, forbidden by the multiperiphcral models, considering 

- Hard Pomeron model 

production of mini-jets. But we will find using Monte Carlo simulations that the pQCD Pomeron is 
unable to adopt the hard channels of hadron production. Then 

- The deep inelastic processes 

for VHM region will be considered to generalize the DGLAP kinematics in the case of heavy QCD jets 
production. The analysis shows that transition to the VHM leads to the necessity to include low- a; 
sub-processes. As a result we get out of the range of pQCD validity. 

- Multiple reduction of jets. 

We will see, that at very high energies in the VHM region the heavy jets creation should be a dominant 
process if the vacuum is stable with reference to the particle production. 

We will consider also decay of the 'false vacuum' to describe the consequence of 

- Phase transition 

in the VHM domain. This channel is hardly seen for n ~ n(s) since the confinement constraints may 
prevent cooling of the system up to phase transitions condition. 



3.1 Peripheral interaction 

3.1.1 Multiperipheral phenomenology 

Later on multiple production physics was developed on the basis of experimental observation (iv). 
The Regge pole model naturally explains this experimental data and, at the same time, absorbs all 
experimental information, (i)-(vi) . At the very beginning, adopting the Regge poles notion without 
its microscopical explanation, this description was self-consistent. The efforts to extend the Regge pole 
model to the relativistic hadron reactions was ended by the Reggeon diagram technique, see [plf and 
references cited therein, and it was used later to construct the perturbation theory for o~ n [ |62| . It was 
shown that the multiplicity distribution is wider then the Poissonian because of the multi-Pomeron 
exchanges. 

The leading energy asymptotics Pomeron contribution reflects the created particle kinematics de- 
scribed in Appendix D, where the available kinematical scenario in the frame of pQCD are described. 
So, the longitudinal momentum of produced particles is large and is strictly ordered. At the same time, 
particles transverse momentum is restricted. 



Let us consider the inelasticity coefficient introduced in (1.2) k = 1 — e rnax /E < 1, where e max 
is the energy of the fastest particle in the laboratory frame. Then the strict ordering of particles in 
the Pomeron kinematics means that n is independent of the index of the particle. So, if the fastest 
particle has the energy e max — (1 — K)y/s, then the following particle should have the energy e\ ~ 
(1 — n)e m ax — (1 — ft) 2 -y/s, and so on. Following to this law, [n — l)-th particle would have the energy 
e„ ~ (1 — K) n y/s. In the laboratory frame the energy should degrade to e„ ~ m. Inserting here the 
above formulated estimation of e n we can find that if the number of produced particles is 

n(s) ~ n £, n = - ln(l - nf > 0, £ = ln(s/m 2 ), (3.1) 

then we may expect the total degradation of energy. This degradation is the necessary condition noting 
that the total cross section of slowly moving particles may depend only slightly on energy and is seems 
necessary for natural explanation of the weak dependence of the hadron cross sections on the energy. 
This consideration would be Lorentz-covariant if one can find the slowly moving particle in an arbitrary 
frame. Resulting estimation of mean mult ipli city have good qualitative experimental confirmation. 

Notice that it was assumed deriving ( |3.l| ) that the energy degrades step by step. In other words, 
if we introduce a time of degradation, then the time ~ £ is needed for complete degradation of energy. 
Assuming the random walk in the normal to incident particle plane, we can conclude that the points 
of particle production are located on a disk (in the moving frame) of radii b ~ £}/ 2 . This means that 
the interaction radii should grow with the energy of the colliding particles. 
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If f(a + b — ► c + ...) is the cross section to observe particle c inclusively in the a and 6 particles 
collision, then it was found experimentally that the ratio 

/(7T+P -> 7T- + ■■■) _ /(/V+p -> 7T~ + ...) _ f{pp-> 7T~ + ...) 

fftot(7T + p) a tot (K+p) <Ttot(pp) 

is universal. This may be interpreted as the direct evidence of fact that the hadron interactions have 
a large-distance character, i.e. that the interaction radii should be large. 

This picture assumes that the probability to have total degradation of energy is 



e 



-6 2 /4a'5 



(3.3) 



where a' is some dimensional constant (the slope of Regge trajectory) and b is the 2-dimensional 
impact parameter. This formulae has also the explanation connected to the vacuum instability with 
reference to the real particle production in the strong colour electric field. 

The above picture has natural restrictions. We can assume that each of the produced particles 
may be the source of above described i-channel cascade of the energy degradation. This means that in 
the frame of the Pomeron phenomenology, we are able to describe the production of 

n < n(s) 2 (3.4) 

particles only. If n > n(s) 2 then the density of particles in the diffraction disk becomes large and (a) 
one should introduce short-distance interactions, or (b) rise interaction radii. It will be shown that 
just (a) is preferable. 

We will build the perturbation theory in the phenomenological frames (i) - (vi). Considering the 
system with variable number of particles the generating function 

T(z, S ) = ]Tz'V n ( S ) (3.5) 

n 

would be useful. One can use also the decomposition: 

T{z, a) = a tot (s) exp Uz- l)d(s) + ~(z - l) 2 C 2 {s) + ... j , (3.6) 



where, by definition, 
is the mean multiplicity, 



C 1 {s)=n{s) = ^-]nT(z,s) 
oz 



C 2 (s)= ^\nT(z,s) 
oz z 



z=l 



(3.7) 
(3.8) 



is the second binomial momentum, and so on. 

Our idea is to assume that all C m , m > 1 may be calculated perturbatively choosing 

P(0,s) = e^ z - 1 ^ s) (3.9) 

as the Born approximation 'superpropagator'. It is evident that ( fj.9| ) leads to Poisson distribution. 
Then, having in mind (ii), (iii) and (v) we will use following anzats: 

P(q, s) = e ai - o) ' a ' q2 ^Se^-Vt, (3.10) 

where the transverse momentum q is conjugate to the impact parameter b. So, the Born term ( |3.1C| ) 
is a Fourier transform of the simple product of (3.9) and (3.3). It contains only one free parameter, 
the Pomeron intercept a(0). On the phenomenological level it is not important to know the dynamical 
(microscopical) origin of ( 3 . 1 0| ) . 
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For our purpose the Laplace transform of P would be useful. If n(s) = ng In s, then 



V(oj,q 2 



r°° 1 

/ <%e-»tp(q 2 ,s)= , (3.11) 
Jo u + a'q 2 + Tp (z) 



It is the propagator of two-dimensional field theory with mass squared 

tpo(z) = (1 — a) + (1 — z)n , n > 0. 

Knowing the Gribov's Reggeon calculus completed by the Abramovski-Gribov-Kancheli (AGK) cutting 
rules (6^] one can investigate the consequences of this approach. 
The LLA approximation of the pQCD Jl^] gives 

A = a(0) - 1 = ^^a s « 0.55, a s = 0.2. (3.12) 

But radiative corrections give A sa 0.2 [Q. We will call this solution as the BFKL model. 

The quantitative origin of the restriction (3.4) is following. The contribution of the diagram with v 
Pomcron exchange gives, since the diffraction radii increase with s, see (3.7), mean value of the impact 
parameter decreasing with v. 

b 2 ~ 4a \n(s/m 2 )/i> — aa— — -, 

where a — A/uq. On the other hand, the number of necessary Pomeron exchanges v ~ n/n(s) since 
one Pomeron gives maximal contribution (with factorial accuracy) at n ~ n(s). In result, 

b "2 ~ aa'^^. (3.13) 

n 

Therefore, if the transverse momentum of created particles is a restricted quantity, i.e. /^gb 2 ~ 1, 
where [io is a constant, then the mechanism of particle production is valid up to 

n~n(s) 2 . (3.14) 

Following our general idea, it will be enough for us find the position of singularity over z. Analysis 
shows that ( 3.12j ) predict the singularity at infinity. 

In Appendix E the Gribov's Reggeon diagram technique with cut Pomerons is described and (3.14) 
is derived. It can be shown using this technique that the model with the critical Pomeron, A = 0, is 
inconsistent from the physical point of view pig ] . 

As was mentioned above, the model with A = a(0) — 1 > is natural for the pQCD. The concrete 
value of A will not be important for us. We will assume only that 

0<A«1. (3.15) 



It is evident that the Born approximation (3.10) with A > violate the Froissart boundary condition. 
But it can shown that the sum of 'eikonal' diagrams^] solves this problem, see [|6| and references cited 
therein. 

The interaction radii may increase with increasing number of produced particles if A > and then 
the restriction (3.14) is not important. In the used eikonal approximation, see Appendix F, 

In 

z ~ z c = 1 + — - In— - (3.16) 

n(s) n(s) 

are essential. Then the interaction radii b 2 ~ B 2 ~ 4c/£(A£ + ln(n/n(s))) for this values of z. Notice 
that 

B 2 - £ 2 >> £ (3.17) 



9 The eikonal approximation in a quantum field theory was developed in 



24 



even for n ~ n max ~ y^i. 

Nevertheless, using (3.16) one can find that the cross section decrease faster than any power of 



- In (^\) = ^ In ^(1 + 0(n( S )/n)) (3.18) 
\crtot(s)J n(s) n(s) 

Generally speaking, although this estimation is right in the VHM region, there may be large corrections 
because of Pomeron self-interactions. But careful analyses shows |63] that these contributions can not 



change drastically the estimation (3.18) 



3.1.2 Dual resonance model 

The search of dynamical source of the Regge description shows the different dynamical nature of 
the Regge and Pomeron poles. The established resonance-Regge pole duality, e.g. |67]] led to the 
Veneziano representation of the Regge amplitudes ||68[ . The Reggeon pole gives the decreasing ~ s~ 1//2 
contribution, but careful investigation shows that the mass spectrum of dual to Regge pole resonances 
increase exponentially. This prediction was confirmed by experiment, see the discussion of this question 

The field theory development is marked by considerable efforts to avoid the problem of colour 
charge confinement. Notice that the classical string has the same excitation spectrum. The remarkable 
attempt in this direction based on the string model, in its various realizations, see e.g. [^0| . But, in spite 
of remarkable success (in formalism especially) there is not an experimentally measurable prediction 
of this approach till now, e.g. [ ffl| . 

We would like to describe in this section production of 'stable' hadrons through decay of resonances 
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Our consideration will use the following assumptions. 

A. The string interpretation of the dual-resonance model bring to the observation that the mass 
spectrum of resonances, i.e. the total number p(m) of mass m resonance excitations, grows exponen- 
tially: 

p(m) = (m/mo) 7 e /3 ° m , 0o = const, m > m . (3.19) 



Note also that the same hadron mass spectrum ( 3. 1S| ) was predicted in the 'bootstrap' approach |39| [74| . 
Moreover, it predicts that 

7 = -5/2. (3.20) 

B. The mass m resonance creation cross section o~ R (m) has the Regge pole asymptotics: 

a R {m) = g R — , g R = const. (3.21) 
m 

It was assumed here that the intercept of the Regge pole trajectory a R = 1/2. So, only the meson 
resonances would be taken into account. 

C. If o~ R (m) describes the decay of a mass m resonance into the n hadrons, then the mean multi- 
plicity of hadrons 

n R {m) = ( 3 - 22 ) 
a R (m) 



Following the Regge model, 



2 

n R (m) = n R \n%. (3.23) 
mi 



H) 

D. We will assume that there is a definite vicinity of n R {m) where a R (m) is defined by n R (m) 
only. So, in this vicinity 

a R {m) = a R {m)e- nIl{ - m) {n R {m)) n /n\. (3.24) 

This is the direct consequence of the Regge pole model, if m/rriQ is high enough. 

Following our idea, we will distinguish the 'short-range' correlations among hadrons and the 'long- 
range' correlations among resonances. The 'connected groups' would be described by resonances and 
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the interactions among them should be described introducing for this purpose the correlation functions 
among strings. So, we will consider the 'two-level' model of hadron creation: the first level describes 
the short-range correlation among hadrons and the second level is connected to the correlations among 
resonances. 

The exact calculations are given in Appendix G. 

Comparing A and B solutions we can see the change of attraction points with rising n: at n ~ 
n 2 (s) = Uq" \n{yfs I mo) the transition from (A) asymptotics to (C) in ( 2.41| ) should be seen. At the 
same time one should see the strong KNO scaling violation at the tail of the multiplicity distribution. 

We have neglected the resonance interactions deriving these results. This assumption seems natural 
since at n(s) « n < h 2 (s) the inequality (G.17) should be satisfied, see discussion of inequalities ( |2. 67 ) 
and (Eggh in Scc.2.3.4. 



3.2 Hard processes 

3.2.1 Deep inelastic processes 

The role of soft colour partons in the high energy hadron interactions is the most intriguing modern 
problem of particle physics. So, the collective phenomena and symmetry breaking in the non-Abelian 
gauge theories, confinement of coloured charges and the infrared divergences of the pQCD are the 
phenomena just of the soft colour particles domain. 

It seems natural that the very high multiplicity (VHM) hadron interaction, where the energy of 
the created particles is small, should be sensitive to the soft colour particle densities. Indeed, the aim 
of this Section is to show that even in the hard by definition deep inelastic scattering (DIS) , see also 
PH | , the soft color particles role becomes important in the VHM region |75| . 

To describe the hadron production in pQCD terms the parton-hadron duality is assumed. This 
is natural just for the VHM process kinematics: because of the energy-momentum conservation law, 
produced (final-state) partons cannot have high relative momentum and, if they were created at small 
distances, production of qq pairs from the vacuum will be negligible (or did not play an important role). 
Therefore, if the 'vacuum' channel is negligible, only the pQCD contributions should be considered 
pc| . f76[ . All this means that the multiplicity, momentum etc. distributions of hadron and colored 
partons are the same. (This reduces the problem practically to the level of QED.) xx 

Let us consider now n particles (gluons) creation in the DIS [[77]. We would like to calculate 
D ab (x, q 2 ;n), where 

53 D ab {x, q 2 ;n) = D ab (x, q 2 ). (3.25) 

n 

As usual, let D ab (x,q 2 ) be the probability to find parton b with virtuality q 2 < in the parton a of 
~ A virtuality, A >> A and a s (X) << 1. We may always choose q 2 and x so that the leading logarithm 
approximation (LLA) will be acceptable. One should assume also that (1/x) >> 1 to have the phase 
space, into which the particles are produced, sufficiently large. 

Then D ab (x, q 2 ) is described by ladder diagrams. From a qualitative point of view this means the 
approximation of random walk over coordinate ln(l/a;) and the time is In In |q 2 |. LLA means that the 
'mobility' ~ ln(l/x)/ lnln \q 2 \ should be large 

hx(l/x) » lnln|g 2 /A 2 | . (3.26) 

But, on other hand [f78[ , 

ln(l/x) << ln|g 2 /A 2 | . (3.27) 

See also Appendix D. 

The leading contributions, able to compensate the smallness of 

a s (A) « 1, 

give integration over a wide range A 2 << k 2 << —q 2 , where k 2 > is the 'mass' of a real, i.e. time-like, 
gluon. If the time needed to capture the parton into the hadron is ~ (1/A) then the gluon should decay 
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if kf >> A 2 . This leads to the creation of (mini)jets. The mean multiplicity nj in the QCD jets is high 
if the gluon 'mass' |fc| is high: lnfL,- ~ ^/ln(k 2 /A 2 ). 

Raising the multiplicity may (i) raise the number of (mini)jets v and/or (ii) raise the mean value 
mass of (mini)jets We will see that the mechanism (ii) would be favorable. 

But raising the me an va lue of gluon masses, \ki\, decreases the range of integrability over fcj, i.e. 
violates the condition Q3.26 ) for fixed x. One can remain the LLA taking x — > 0. But this may 
contradict to (3.27), i.e. in any case the LLA becomes invalid in the VHM domain and the next to 
leading order corrections should be taken into account. 

Noting that the LLA gives the main contribution, that the rising multiplicity leads to the infrared 
domain, where the soft gluon creation becomes dominant. 

First of all, neglecting the vacuum effects, we introduce definite uncertainty to the formalism. It is 
reasonable to define the level of strictness of our computations. Let us introduce for this purpose the 
generating function T ab {x, q 2 ; z): 

1 f dz 

D ab (x, q 2 ;n) = — j> -^T ab {x, q 2 ; z). (3.28) 

At large n, the integral may be calculated by the saddle point method. The smallest solution z c of the 
equation 

n = z-^-bxT ab (x,q 2 ;z) (3.29) 
az 

defines the asymptotic over n behavior: 

D ab (x, q 2 ;n) ex exp{— rt \nz c (x, q 2 ; n)}. (3.30) 

Using the statistical interpretation of z c as the fugacity it is natural to write: 

2 C ab (x,q 2 ;n) 
mz c (x, q ; n) = — — ^— . (3.31) 

n ab (x,q 2 ) 

Notice that the solution of eq.( 3.2S| ) z c (x, q 2 ; n) should be an increasing function of n. At first glance 
this follows from the positivity of all D ab (x,q 2 ;n). But actually this assumes that T ab (x,q 2 ; z) is a 
regular function of z at z = 1. This is a natural assumption considering just the pQCD predictions. 
Therefore, 

D ab (x, q 2 ;n) cx exp{-^ — ^-C ab (x, q 2 ;n)}. (3.32) 

n ab (x,q 2 ) 

This form of D ab (x, q 2 ; n) is useful since usually C ab (x, q 2 ; n) is a slowly varying function of n. So, for 
a Poisson distribution C ab (x, q 2 ;n) ~ Inn. For KNO scaling we have C ab (x, q 2 ;n) — const, over n. 

We would like to note that, neglecting effects of vacuum polarization, we introduce into the ex- 
ponent so high uncertainty assuming n ~ n p that it is reasonable to perform the calculations with 
exponential accuracy. So, we would calculate 

- p,ab{x,q 2 ;n) = In D °*( x ' q L^j = _ — " — C ab (x, q 2 ; n)(l + 0(l/n)) (3.33) 
D ab (x,q 2 ) n ab (x,q 2 ) 

The n dependence of C ab (x,q 2 ;n) defines the asymptotic behavior of p, ab (x,q 2 ;n) and calculation 
of its explicit form would be our aim. 

We can conclude, see Appendix H, that our LLA is applicable in the VHM domain till 

uj(t, z) « ln(l/a:) << r = ln(-g 2 /A), (3.34) 

where 

oj(t,z)= f '*L w s( T >, 2 ). (3.35) 



and 



w a (r,z) = Vz"<(r) (3.36) 
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is the generating function of the multiplicity distribution in a gluon jet. In the frame of constraints 

fH), 

F ab (q 2 ,x;w) oc exp ^Nuj{t, z) ln(l/x)| . (3.37) 



The mean multiplicity of gluons created in the DIS kinematics 
d 



n g {r,x) = ^-\nF ab (q 2 ,x;w) |, =1 = AN /x)/lnr » wi(r), (3.38) 
where 

Wl (r)= f — %(r) (3.39) 
and the mean gluon multiplicity in the jet %(t) has the following estimation J79|: 

lnnj(r)~Vr (3.40) 

Inserting ( |3.40| ) into ( |3.39| ), 



Wl (r) = nj(r)/Vf. 



Therefore, noting (3.27) 



%(t, a;) ~ rij (^^^^(l/^/rlnr << %(r). (3.41) 
This means that the considered 't-channel' ladder is important in the narrow domain of multiplicities 

n~ n g «n 3 . (3.42) 

So, in the VHM domain n >> n g one should consider 

(i) The ladder diagrams with a small number of rungs; 

(ii) To take into account the multi-jet correlations assuming that increasing multiplicity leads to the 
increasing number of rungs in the ladder diagram. 

To choose one of these possibilities one should consider the structure of u>(t, z) much more carefully. 
This will be done in the Sec. 3. 5. 

We can conclude that in the VHM domain, multiplicity production unavoidably destroys the ladder 
LLA. To conserve this leading approximation one should choose x —> and, in result, to get to the 
multi-ladder diagrams, since in this case a s ln(— q 2 /\ 2 ) ~ 1 and a s ln(l/x) ~ 1. Such theory was 
considered in |p0| . 

3.2.2 QCD jets 

As was mentioned above, the pQCD description is right if the colour particles virtuality is bounded 
from below, \q 2 \ > la 2 , where A is chosen so that a s {\ 2 ) « 1. This kinematical restriction leads to 
the infrared cutoff |8l], ^2| and may essentially influence the particle production in the VHM region. 
It is a special property of pQCD. Indeed, for example, careful investigation of this question in the 
asymptotically free (c/? 3 )6-theory ||| shows that this restriction is 'unobservable' since their inclusion 
takes us beyond the LLA (^4). At the same time, the condition \q 2 \ > la 2 essentially shrinks the phase 
space where particles are produced. 

Particle (gluons) distribution in pQCD jets was investigated firstly in |^l[ |8f| and it was shown 
that the generating function is singular at z s — 1 ^ (1/n). Let us consider this solution stable with 
reference to the discussed cutoff. 

The explicit formulae for one jet production may be written in the form, see Appendix I: 

o£ y) (M) = a^(M, n)e- c i n ' n ^ M \ n > %(M), (3.43) 

where a^\n,M) is the polynomial function of n, fij(M) is the mean multiplicity in the mass M jet 
and Cj is a positive constant. 
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The linear behavior of the exponent in {1.4) over n/n has important consequences. So, let us assume 
that the total energy M is divided into two jets of masses M\ and Mi equally: Mi = Mi — M/2. 
If, for instance, Mi << M\ ~ M then the distribution will coincide with (3.43), but the second jet 
distribution would renormalize the coefficient o/ 1 ^ . 

Then the multiplicity distribution in the two-jet event would be 



a^\M) = a^ j \M,n)e- c ^ n ^ 2 \ 

where n\ + m = n is the total multiplicity. 

Comparing (3.43) with (3.44) we can see that with exponential accuracy: 



(3.44) 



exp{- 



n 3 {M)-n,{M/2) 
' fi j {M)n j {M/2) 



the (HH ) would dominate in the VHM domain since the mean multiplicity rij{M) increases with M. 

The experimental observation of this phenomena crucially depends on the value of a J , a 2 -?,... but 
if ( 3.43|) is satisfied then one can expect that the events in the VHM domain would be enhanced by 
QCD jets and the mass of jets would have a tendency to be high with growing multiplicity. 

The singular at finite z solutions arise in the field theory, when the s-channel cascades (jets) are 
described ]5q] . By definition T{z,s) coincides with the total cross section at z = 1. Therefore, the 
nearness of z c to one defines the significance of the corresponding processes. It is evident that both s 
and n should be high enough to expect the jets creation. 

Summarizing the above estimations, we may conclude that 



0{e~ n ) < a n < 0(l/n), 



(3.45) 



i.e. the soft Regge-like channel of hadron creation is suppressed in the VHM region in the high energy 
events with exponential accuracy. 



3.3 Phase transition - condensation 

The aim of this section is to find the experimentally observable consequences of collective phenomena in 
the high energy hadron inelastic collision fl86| . We will pay the main attention to the phase transitions, 
leaving out other possible interesting collective phenomena. 

The statistics experience dictates that we should prepare the system for the phase transition. The 
temperature in a critical domain and the equilibrium media are just these conditions. It is evident that 
they are not a trivial requirement considering the hadrons inelastic collision at high energies. 

The collective phenomena by definition suppose that the kinetic energy of particles of media are 
comparable, or even smaller, than the potential energy of their interaction. It is a quite natural condi- 
tion noting that, for instance, the kinetic motion may destroy even completely at a given temperature 
T, necessary for the phase transition long-range order. This gives, more or less definitely, the critical 
domain. 

The same idea as in statistics seems natural in the multiple production physics. We will assume 
(A) that the collective phenomena should be seen just in the very high multiplicity (VHM) events, 
where, because of the energy-momentum conservation laws, the kinetic energy of the created particles 
can not be high. 

We will lean at this point on the S'-matrix interpretation of statistics |]87| , see Sec. 2. 2. It based on 
the S'-matrix generalization of the Wigner function formalism of Carruthcrs and Zachariazen p9| and 
the real-time finite temperature field theory of Schwinger and Keldysh |2^, ^30|, see Appendices A,B 
and C. 

It was mentioned that the n-particle partition function in this approach coincides with the n particle 
production cross section a n {s) (in the appropriate normalization condition). Then, the cross section 
<T n {s) can be calculated applying the n-point Wigner function W n {Xi, Xi, ...,X n ). In the relativistic 
case Xk = {u, q)k are the 4- vectors. So, the external particles are considered as the 'probes' to measure 
the state of the interacting fields, i.e. the low mean energy of probes means that the system is 'cold'. 
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The multiple production phenomena may be considered also as the thermalization process of in- 
cident particle kinetic energy dissipation into the created particle mass. From this point of view the 
VHM processes are highly noncquilibrium since the final state of this case is very far from the initial 
one. It is known in statistics [^4j that such processes aspire to be the stationary Markovian with a high 
level of entropy production. In the case of complete thermalization, the final state is in equilibrium. 

The equilibrium we will classify as the condition in the frame of which the fluctuations of corre- 
sponding parameter are Gaussian. So, in the case of complete thermalization, the probes should have 
the Gaussian energy spectra. In other terms, the necessary and sufficient condition of the equilibrium 
is the smallness of the mean value of energy correlators |87j . From the physical point of view, the 
absence of these correlators means depression of the macroscopic energy flows in the system. 

The multiple production experiment shows that the created particle energy spectrum is far from a 
Gaussian law, i.e. the final states are far from equilibrium. The natural explanation of this phenomena 
consist in the presence of (hidden) conservation laws in the interacting Yang-Mills fields: it is known 
that the presence of sufficient number of first integrals in involution prevents thermalization completely. 

But nevertheless the VHM final state may be equilibrium (B) in the above formulated sense. This 
means that the forces created by the non-Abelian symmetry conservation laws may be frozen during 
the thermalization process (remembering its stationary Markovian character in the VHM domain) . We 
would like to take into account that the entropy S of a system is proportional to number of created 
particles and, therefore, S should tend to its maximum in the VHM region 0. 

One may consider following the small parameter (n(s)/n) << 1, where n(s) is the mean value of 
the multiplicity n at a given CM energy ^fs. Another small parameter is the energy of the fastest 
hadron e max . One should assume that in the VHM region (e m ax/\^s) — > 0. So, the conditions: 

n(£) << ernax ^ Q . 
71 V s 

would be considered as the mark of the processes under consideration. We can hope to organize the 
perturbation theory over them having there small parameters. In this sense VHM processes may be 
'simple', i.e. one can use for their description semiclassical methods. 

So, considering VHM events one may assume that the conditions (A) and (B) are satisfied and 
one may expect the phase transition phenomena. 

The S-matrix interpretation of statistics is based on the following definitions. First of all, let us 
introduce the generating function Q : 

T(z,8)=^z n a n {s). (3.47) 

n 

Summation is performed over all n up to n max — -y/s/m and, at finite CM energy y/s, T(z,s) is a 
polynomial function of z. Following our idea, see Sec. 2. 3, let us assume now that z > 1 is sufficiently 
small and for this reason T(z, s) depends on the upper boundary n max only weakly. In this case one 
may formally extend summation up to infinity and in this case T(z, s) may be considered as a whole 
function. This possibility is important being the equivalent of the thermodynamical limit and it allows 
to classify the asymptotics over n in accordance with the position of singularities over z. 

Let us consider T(z) as the big partition function, where z is 'activity'. It is known |3lf that T(z) 
should be regular inside the circle of unit radius. The leftist singularity lies at z = 1. This singularity 
is manifestation of the first order phase transition |5l|, |24|, |5^] . 

The origin of this singularity was investigated carefully in the paper |p4j| . It was shown that the 
position of singularities over z depends on the number of particles n in the system: the two complex 
conjugated singularities move to the real z axis with rising n and in the thermodynamical limit, n — oo, 
they pinch the point z = 1 in the first order phase transition case. More general analysis [[32] shows 
that if the system is in equilibrium, then T(z) may be singular only at z = 1 and z = oo. 

The position of the singularity over z and the asymptotic behavior of a n are related closely. Indeed, 
for instance, inserting into ( |3.47 ) a n oc exp{— en 1 } we find that T(z) is singular at z = 1 if 7 < 1. 



Generally, using the Mellin transformation (2.37) one can find an asymptotic estimation ( 2.39| ) 



-n In z c {n) 



z c > 1, (3.48) 
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where z c is smallest solution of the equation of state 



d 

n = z—\nT{z). (3.49) 

Therefore, to have the singularity at z = 1, we should consider z c (n) as a decreasing function of n. 
On the other hand, at constant temperature, lnz c (n) ~ Hc{n) is the chemical potential, i.e. is the work 
necessary for creation of one particle. So, the singularity at z = 1 means that the system is unstable: 
the less work is necessary for creation of one more particle if /J.(n) is a decreasing function of n. 



The physical explanation of this phenomena is the following, see also 1 88 . The generating function 
T(z) has following expansion: 

T{z) =exp{^A}, (3.50) 

where bi are known as the Mayer's group coefficients They can be expressed through the inclusive 
correlation functions and may be used to describe formation of droplets of correlated particles, see 
Sec. 2. 3. 3. So, if droplet consist of I particles, then 

k ~ e -W (d - 1)/d (3.51) 

is the mean number of such droplets. Here is the surface energy of d-dimensional droplet. 

Inserting this estimation into ( 3.50| ), 



lnT(z)~^ e ^-«^ 1,/d ), ^ = lnz. (3.52) 
l 



The first term in the exponent (31 fi is the volume energy of the droplet and being positive it tries to 
enlarge the droplet. The second surface term — tries to shrink it. Therefore, the singularity 
at z = 1 is the consequence of instability: at z > 1 the volume energy abundance leads to unlimited 
growth of the droplet. 

In conclusion we wish to formulate once more the main assumptions. 

(I) . It was assumed first of all that the system under consideration is in equilibrium. This condition 
may be naturally reached in the statistics, where one can wait the arbitrary time till the system 
becomes in equilibrium. Note, in the critical domain, the time of relaxation t r ~ (T c / (T — T C )Y — > oo, 
(T — T c ) — > +0, v > 0, T c is the critical temperature. 

We cannot give the guarantee that in the high energy hadron collisions the final state system is 
in equilibrium. The reason of this uncertainty is the finite time the inelastic processes and presence 
of hidden (confinement) constraints on the dynamics. But if the confinement forces are frozen in the 
VHM domain, i.e. the production process is 'fast', then the equilibrium may be reached. 

We may formulate the quantitative conditions, when the equilibrium is satisfied [42]] . One should 
have the Gaussian energy spectra of created particles. If this condition is hardly investigated in the 
experiment, then one should consider the relaxation of 'long-range' correlations. This excludes the 
usage of relaxation condition for the 'short-range' (i.e. resonance) correlation 

(II) . The second condition consists in the requirement that the system should be in the critical 
domain, where the (equilibrium) fluctuations of the system become high. Having no theory of hadron 
interaction at high energies we can not define where lies the 'critical domain' and even whether it exists 
or not. 

But, having the VHM 'cold' final state, we can hope that the critical domain is achieved. 
The quantitative realization of this picture is given in Appendix J. It is important to note that 
used there the semiclassical approximation is rightfully in the VHM domain. 



4 Conclusion 
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4.1 Discussion of physical problems 



It seems useful to start the discussion of models outlined the main problems, from authors point of 
view. 

A. Soft colour parton problems 

The infrared region of soft colour parton interactions is a very important problem of high energy 
hadron dynamics. Such fundamental questions as the infrared divergences of pQCD, collective phe- 
nomena in the coloured particles system and symmetry breaking are the phenomena of the infrared 
domain. 

The standard (most popular) hadron theory considers pQCD at small distances (in the scale of 
A ~ 0.2Gev) as the exact theory. This statement is confirmed by a number of experiments, namely 
deep-inelastic scattering data, hard jets observation. But the pQCD predictions have a finite range of 
validity since the non-perturbative effects should be taken into account at distances larger then 1/A. 

It is natural to assume, building the complete theory, that at large distances the non-perturbative 
effects lay or^\ the perturbative ones. As a result pQCD loses its predictability screened by the 
non-perturbative effects. 

Notice, the pQCD running coupling constant a s (q 2 ) = l/61n(g 2 /A 2 ) becomes infinite at q 2 = A 2 
and we do not know what happens with pQCD if q 2 < A 2 . There are few possibilities. For instance, 
there is a suspicion jl3j that at q 2 ~ A 2 the properties of theory changed so drastically (being defined 
on a new vacuum) that even the notions of pQCD disappeared. This means that pQCD should be 
truncated from below on the 'fundamental' scale A. It seems natural that this infrared cut-off would 
influence the soft hadrons emission. 

The new possibility is described in Appendix K. This strict formalism allows to conclude that pure 
pQCD contributions are realized on zero measure, i.e. it is the phenomenological theory only. The 
successive approach shows that the Yang-mills theory should be described in terms of (action, angle)-like 
variables. The last one means that the self-consistent description excludes such notions as the 'gluon'. 
As a result of this substitution new perturbation theory would be free from infrared divergences, i.e. 
there is not necessity to introduce the infrared cutoff parameter A. (Moreover, in the sector of vector 
fields (without quarks) the theory is ultraviolet stable.) It seems important for this reason to investigate 
experimentally just VHM events, where the soft colour partons production is dominant. 

It is important to try to raise the role of pQCD in the 'forbidden' area of large distances. The 
VHM processes are at highly unusual condition, where the non-perturbative effects must be negligible. 

B. Dissipation problems 

The highly nonequilibrium states decay (thermalization) which means in the pQCD terms that the 
process of VHM formation should be enhanced, at least in asymptotics over multiplicity and energy, 
by jets. It is the general conclusion of nonequilibrium thermodynamics and it means that the very 
nonequilibrium initial state tends to equilibrium (thermalized) as fast as possible. 

The entropy S of a system is proportional to the number of created particles and, therefore, iS 
should tend to its maximum in the VHM region H. But the maximum of entropy testifies also the 
equilibrium of the system. 

C. Collective phenomena 

We should underline that the collective phenomena may take place if and only if the particles 
interaction energy is comparable to the kinetic one. The VHM system considered may be 'cold' 
and 'equilibrium'. For this reason the VHM state is mostly adopted for investigation of collective 
phenomena. One of possible states in which the collective effects, see a.g. [^9|, may be important is 
the 'could coloured plasma' pQ |. 

The fundamental interest presents the problem of vacuum structure of Yang-Mills theory. For 
instance, if the process of cooling is 'fast', since the dissipation process of VHM final state forma- 
tion should be as fast as possible, then one may consider formation of vacuum domains with various 
properties. Then decay of these domains may lead to large fluctuations, for instance, of the isotopic 
spin. 

10 The corresponding formalism was described e.g. in Jig}] , 
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Another important question is the collective phenomena in the VHM final state. The last one 
may be created 'perturbatively', for instance, by the formation of heavy jets. Then the colour charges 
should be confined. There are various predictions about this process. One of them predicts that there 
should be a first order phase transition. 

4.2 Model predictions 

Now we can ask: what can models say concerning the above problems? 

A. Soft coloured partons production 

The multiperipheral models predict fast decreasing of topological cross section in the VHM domain 
n(s) 2 << n << n max , a n < 0(e~ n ). At the same time the mean transverse momentum should 
decrease in this domain since the interaction radii should 'increase' with n. 

The BFKL Pomeron predicts the same asymptotics, a n < 0(e~ n ), but the pQCD jet predicts 
o~ n = 0(e~ n ). The naive attempt to insert into the BFKL Pomeron the production of particles via 
(mini)jets seems impossible. 

This 'insertion' can be done into the DIS ladder but investigation of the LLA kinematics in the 
VHM domain allows the conclusion that the 'low-x' contributions should be taken into account. 

All this experience allows the assumption that in the VHM domain no 'i-channel ladder' diagrams 
play sufficient role. This means the existence of a transition to the processes with jet dominance. The 
pQCD is unable to predict the transition mechanism. 

B. Transition into 'equilibrium' 

If the t-channel ladders are 'destroyed' in the VHM region, then jets, despite the small factor 0(1/ s) 
in the cross sections, are the only mechanism of particle production in the VHM domain. Dominance 
of heavy jets in the VHM domain may naturally explain the tendency towards equilibrium. 

But the description of thermalization in terms of jets of massless gluons production destroy this 
hope: the jet contribution cr„ = 0(e~ n ) assumes 'bremsstrahlung' of soft gluons 0. This prevents 
the equilibrium since ordering without fail introduce the non-relaxing correlations. 

C. Collective phenomena 

Considering the collective phenomena, we proposer to distinguish (a) the collective phenomena 
connected with the vacuum and (b) the collective phenomena produced in the VHM system. Following 
the experience of Sec. 3. 3 we can conclude that the signal of vacuum instability is inequality: a n > 
0{e~ n ). 

The case (b) will not effect the cross section a n . But if the system became equilibrium in the VHM 
domain then the collective phenomena may be investigated using ordinary thermodynamical methods. 
For instance, noting that — {ln(T(/3 c , z)/T(/3 c , z))}/(5 c = T((i c ,z) is the free energy one can measure 
the thermal capacity 

^-F(f3 c ,z) = C(f3 c ,z). (4.1) 

Then, comparing capacities of hadron and 7-quanta systems we can say whether or not the phase 
transition happened. 

The connection of the equilibrium and relaxation of correlations is well known ^2j. Continuing 
this idea, if the VHM system is in equilibrium one may assume that the colour charges in the pre- 
confinement phase of VHM event form the plasma. One should note here that the expected plasma is 
'cold' and 'dense'. For this reason no long-range confinement forces would act among colour charges. 
Then, being 'cold', in such a system various, collective phenomena may be important. 

4.3 Experimental perspectives 

The experimental possibilities in the VHM domain are not clear up to now. Nevertheless first steps 
toward formulation of trigger system was done, see |)2j. By this reason we would like to restrict 
ourselves by following two general questions. It seems that these questions are mostly important being 
in the very beginning of VHM theory. 
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I. For what values of multiplicity at a given energy the VHM processes become hard! 

The answer to this question depends on the value of the incident energy. If we know the answer then 
it will appear possible to estimate 

- the role of multiperipheral contributions, 

- the jet production rate, 

- the role of vacuum instabilities. 

It seems that the experimental answer to this question is absent since produced particles are soft, 
theirs mean transverse and longitudinal momenta have the same value. In our understanding this 
question means: the total transverse energy may be extremely high. 

It is interesting also to search the heavy jets, i.e. to observe the fluctuations of particle density in 
the event-by-event experiment, but this program seems vague since, for all evidence, fractal dimensions 
tend to zero with increasing multiplicities. 

II. For what values of multiplicity does the VHM final state become equilibrium? 
We hope that having answer on this question we would be able 

- to investigate the status of pQCD, 

- to observe the phase transition phenomena directly, 

- to estimate the role of confinement constraints. 

The equilibrium means that the energy correlation functions mean values are small. It is interesting 
also, for example, the charge equilibrium, when the mean value of charge correlation functions are small. 

Notice, the effect of the phase space boundary may lead to 'equilibrium'. Indeed, if e(p) ~ m + 
p 2 /2m and p 2 << m 2 then one may neglect the momentum dependence of the amplitudes a n . In this 
case the momentum dependence is defined by the Boltzmann exponent e _/3p2 / 2m , {3 — > oo, only and 
we get naturally to the Gaussian law for momentum distribution. Correlators should be small in this 
case since there is no interactions among particles (a n are constants). But our question assumes that 
we investigate the possibility of equilibrium when n << n max and p 2 >> m 2 . 



A Appendix. Matsubara formalism and the KMS 
boundary condition 

There are various approaches to build the real-time finite-temperature field theories of Schwinger- 
Keldysh type (e.g. Ej). All of them use various tricks for analytical continuation of imaginary-time 
Matsubara formalism to real time |)3j . The basis of the approaches is the introduction of the Matsubara 
field operator 

$ M (x,/3)=e^<fs(x) e -^, (A.l) 
where ( l > s(x) is the interaction-picture operator introduced instead of the habitual Heisenberg operator 

*(x,t) = e ltH <$>s{*)e- ltH . 



Eq.( A.l ) introduces the averaging over the Gibbs ensemble instead of averaging over zero-temperature 
vacuum states. 

If the interaction switched on at the moment tj adiabatically and switched off at tf then there is 
the unitary transformation: 

$(x, t) = U(ti,t f )U(ti,t)$s(x)U(t, U). (A.2) 

Introducing the complex Mills time contours (4(| to connect ti to t, t to tf and tf to U we form a 
'closed-time' contour C (the end-points of the contours joined together). This allows to write the last 



equality (A.2) in the compact form: 

*(*)= T c {$(x)e l /c dVi ""^')} s , 
where Tc is the time-ordering on the contour C operator. 
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The generating functional Z(j) of correlation (Green) functions has the form: 



Z(j) = R(0) < Tce '/ ^W*)+j(*)*(*)} S > 

where <> means averaging over the initial state. 

If the initial correlations have a little effect, we can perform averaging over the Gibbs ensemble. 
This is the main assumption of the formalism: the generating functional of the Green functions Z(j) 
has the form in this case: 



Z(j) = /W < <f / ;t l \e-" H T c eIc d4x ^^\^;t l 



> 



with $' = <f>'(x). In accordance with (AT) we have: 

<$';t i |e-^=<$ / ; U - i(3\ 

and, as a result, 

I" if d 4 x{L(x)+j(x)<S>(x)} 

Z(j) = / D&e Jc ? (A.3) 
where the path integration is performed with KMS periodic boundary condition: 

= i/3). 



In ( A.3 ) the contour C@ connects tj to tt, tf to ij and ij to ti — i[3. Therefore it contains an imaginary- 
time Matsubara part ti to t% — i(3. A more symmetrical formulation uses the following realization: ti to 
tf, tf to tf — i(3/2, tf — i/3/2 to t% — ifl/2 and ti — i/3/2 to t j — i /3 (e.g. J|^]). This case also contains the 



imaginary-time parts of the time contour. Therefore, eq. (A^) presents the analytical continuation of 
the Matsubara generating functional to real times. 



One can note that if this analytical continuation is possible for Z(j) then representation ( A.3 ) gives 
good recipe of regularization of frequency integrals in the Matsubara perturbation theory, e.g. |47f . 
But it gives nothing new for our problem since the Matsubara formalism is a formalism for equilibrium 
states only. 

B Appendix. Constant temperature formalism 

The starting point of our calculations is the n- into m-particles transition amplitude a nm , the derivation 
of which is well known procedure in the Lehmann-Symanzik-Zimmcrmann (LSZ) reduction formalism 
[9^ ] framework, see also fl95| . The (n + m)-point Green function G nm are introduced for this purpose 
through the generating functional Zj J9t| : 

n m 

G nm (x, y) = H)" +r " II ft**) II Kvk)Zj, (B.l) 



fc=l fc=l 



where 



and 



The action 



Sj{x) 



Zj = / ZX&e 15 ^*). (B.3) 



S^) - S m ~ V($) + fdxj(x)${x), (B.4) 

where Sq($) is the free part and V($) describes the interactions. At the end of the calculations one 
can put j = 0. 
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To provide the convergence of the integral (B.3) over the scalar field $ the action Sj(§>) must 
contain a positive imaginary part. Usually for this purpose Feynman's «e-prescription is used. But it 
is better for us to use the integral on the Mills complex time contour C+ jl6| ^|. For example, 

C± : t^t + ie,e^+Q, - oo < t < +00 (B.5) 

and after all the calculations return the time contour on the real axis putting e = 0. 

In eq. (B.3) the integration is performed over all field configurations with standard vacuum bound- 
ary condition: 

J d i xd^d^^>) = J dcr^d^® = 0, (B.6) 

which leads to zero contribution from the surface term. 
Let us introduce now field <f> through the equation: 



and perform the shift $ — > $ + cp in integral ( |B.3| ), conserving boundary condition (B.6). Considering 
4> as the probe field created by the source: 



4>{x) = J dyG (x - y)j(y), 
(d 2 + m%G (x-y) =6(x-y), (B. 
only the connected Green function G c nm will be interesting for us. Therefore, 



G c nm {x, y) = H)" +m J] j(x k ) ]J j{Vk)Z{<P), (B.9) 
fc=i fc=i 

where 

Z((f>) = J Z?$e 4S W- iV (*+« (B.10) 

is the new generating functional. 

To calculate the nontrivial elements of the S'-matrix we must put the external particles on the mass 
shell. Formally this procedure means amputation of the external legs of G c nm and further multiplication 
on the free particle wave functions. As a result the amplitude of n- into m-particles transition a nm in 
the momentum representation has the form: 

n m 

a nm (q,p) = H) n+ ™ II II P(PkM4>)- (B.11) 
k=l fc=l 

Here we introduce the particle distraction operator 

<Mg) = f dxe-^^x), fa) = (B.12) 
J 54>(x) 

Supposing that the momentum of particles are insufficient for us the probability of n- into m- 
particlcs transition is defined by the integral: 

r nm = —. — ; / dn n (q)dQ m (p)5 (i) (V \q k - y"p fc )|a„ m | 2 , (B.13) 
n to / * — ' * — ' 



k=l k=l 



where 



k=l fe=l v ; KHKJ 
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is the Lorentz-invariant phase space clement. We assume that the energy- momentum conservation 
5-function was extracted from the amplitude. 

Note that r nm is the divergent quantity. To avoid this problem with trivial divergence, connected 
integration over reference frame, let us divide the energy-momentum fixing ^-function into two parts: 

<5 (4) (E*-Ep*) = /^ (4) (^-E^ (4) ( p -E^) (B.15) 
and consider a new quantity: 

-. « n n 

R ( p ) = y2-n dn n (q)dn m ip)s^Hp-y^q k )s^\p-y^ Pk )\a nm \ 2 (B.ie) 

L — ' n\m\ / £ — ' £ — ' 

n,m J k=l k=l 



defined on the energy momentum shell (2.6). Here we suppose that the numbers of particles are not 
fixed. It is not too hard to see that, up to a phase space volume, 



R = J d A P R{P) (B.17) 



is the imaginary part of the amplitude < vac\vac >. Therefore, computing r(P) the standard renor- 
malization procedure can be applied and the new divergences will not arise in our formalism. 
The Fourier transformation of (5-functions in (B.16) allows one to write R(P) in the form: 

R < p » - / (R18) 



where 



1 r m n 

=Es Y[{dn(q k )e- ia ^} l[{dn(p k )e- ia ^}\a nm \ 2 . (B.19) 

ir* im ' ' " 7 1 7, 1 



fc=l k=l 



Introducing the Fourier-transformed probability p(a) we assume that the phase-space volume is 
not fixed exactly, i.e. it is proposed that the 4-vector P is fixed with some accuracy if cti are fixed. 
The energy and momentum in our approach are still locally conserved quantities since an amplitude 
a nm is translationally invariant. So, we can perform the transformation: 

"i E qk = ( ai _ E 9fe + ai E 9fc ~* ( ai ~ cri ) E qk + aiP (B.20) 

since 4-momenta are conserved. The choice of <j\ fixes the reference frame. This degree of freedom of 
the theory wa s cons idered in [Tl| . 

Inserting flB.lip into (|B.19|) we find that 

p(a) = exp{— i J dxdx' (4>+ (x) (a; — x', q.2)4>— (x') — 

-4-(x)D-+(x - x , ,a 1 )$+(a/))}Z(tf,+)Z*{-4>-), (B.21) 
where and D |_ are the positive and negative frequency correlation functions correspondingly: 

^f-^^/ft^^ (B,2, 

describes the process of particles creation at the time xq and its absorption at x' Q , xq > x' , and a is 
the CM 4-coordinate. Function 

D-+(x-x',a) = i J ' dn 1 (q)e- iq{x - x ' +a) (B.23) 

describes the opposite process, xq < x' . These functions obey the homogeneous equations: 

{d 2 + m 2 ) x G+- = (d 2 + m 2 ) x G-+ = (B.24) 
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since the propagation of on mass-shell particles is described. 

We suppose that Z(<f>) may be computed perturbatively. For this purpose the following transfor- 
mations will be used (X = 8/5X at X = 0): 



-iV(</>) _ g-» J dxj(x)(ji'(x) e i J dxj(x)4>(x) e -iV(4>') 

_ e J dx</,(x)4>'(x) e ^ l v(4,') _ 

_ e -iV(-ij) e i J dxj(x)<p(x) 



(B.25) 



where j was defined in (BJ2) and <f> in ( B.12 ). At the end of the calculations, the auxiliary variables j, 
4>' can be taken equal to zero. Using the first equality in ( |B.25| ) we find that 



Z((j>) = J dx K x )H x ) g-«^(*+0) g-| J dxdx'j(x)D ++ (x-x')j{x') 

where -D++ is the causal Green function: 

(d 2 + m 2 ) x G ++ (x-y) = 6(x-y). 



(B.26) 



(B.27) 



Inserting ( B.26| ) into (B.2I ) after simple manipulations with differential operators, see (B.25), we may 
find the expression: 



p(a) 



-iV(-ij+)+iV(-ij-) 



exp{^ J dxdx' 



where 



x (.U( x ) D +-(x- x',a{)j-{x') - j-(x)D- + (x - x',a 2 )j+(x') 
-j+(x)D ++ (x ~ x')j + {x') + j-(x)D— (x - x')j-(x'))}, 

D— = (D++Y 



(B.28) 
(B.29) 

is the anticausal Green function. 

Considering the system with a large number of particles, we can simplify the calculations choosing 
the CM frame P = (Pq = E,0). It is useful also JO], |3j| to rotate the contours of integration over 

«o,fc : ao,k = -«/3fc,Im/3 fe = 0, k= 1,2. 

For the result, omitting the unnecessary constant, we will consider p — p{0). 

External particles play a double role in the S'-matrix approach: their interactions create and 
annihilate the system under consideration and, on the other hand, they are probes through which the 
measurement of a system is performed. Since /3k are the conjugate to the particles energies quantities 
we will interpret them as the inverse temperatures in the initial and final (/?2) states of interacting 
fields. They are the 'good' parameters if and only if the energy correlations are relaxed. 



Kubo-Martin-Schwinger boundary condition 



The simplest (minimal) choice of Q(<Jco) assumes that the system under consideration is 
surrounded by black-body radiation. This interpretation restores Niemi-Semenoff 's formulation of the 
real-time finite temperature field theory [p8[ . 

Indeed, as follows from (B.21), the generating functional p(ot) is defined by corresponding gener- 
ating functional 



pot 



6±) = Z(ci) + )Z*(-<j)-) = JZ3$ + D$_e iS o('i'+)-^o(*-) x 



xe 



-iV{<t> + +4> + )+iV(<s> --<t>-) 



(B.30) 



see ( B.21 ). The fields (0+,$+) an d (tft_ , $_) were defined on the time contours C+ and C_. 

As was mentioned above, see ( [2.11 ), the path integral ( B.30 ) describes the closed path motion in 
the space of fields <E>. We want to use this fact and introduce a more general boundary condition which 
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also guarantees the cancelation of the surface terms in the perturbation framework. We will introduce 
the equality: 

' dcr At $+d' < $ + = f da^-d**®-. (B.31) 



The solution of eq.(B.31) requires that the fields <&+ and $_ (and their first derivatives d fi &±) coincide 
on the boundary hypersurface a^: 

$±(o-oo) = ®M 7^0, (B.32) 

where, by definition, ^((Too) is the arbitrary 'turning-point" field. 

In the absence of the surface terms, the existence of a nontrivial field $(<7oo) has the influence only 
on the structure of the Green functions 

G++ =< >, G+_ =< >, 

G-+ =< >, G__ =< T$_$_ >, (B.33) 

where T is the antitemporal time ordering operator. These Green functions must obey the equations: 

(8 2 + m 2 ) x G+- (x-y) = (8 2 + m 2 ) x G-+ (x - y) = 0, 
(d 2 + m 2 ) x G++{x -y) = (d 2 + m 2 )*G__(z - y) = S(x - y), (B.34) 

and the general solution of these equations: 

Gii — Da ~~\~ Qa 3 

Gij = 9ij-, i .i (B.35) 
contain the arbitrary terms which are the solutions of homogenous equations: 

(d 2 + m 2 ) x g l3 (x~y) = 0, i,j = +,-. (B.36) 
The general solution of these equations (they are distinguished by the choice of the time contours C±) 

9ij (x - x') = I ' d^iqy^-^mM (B.37) 



are defined through the functions which are the functionals of the 'turning-point' field $(croo): if 
< i ) (o'oo) = we must have Ujj = 0. 

Our aim is to define riij. We can suppose that 

Uij ~< $(ctoo)---*(o-oo) > ■ 

The simplest supposition gives: 

n i:j ~< >~< $ 2 (aoo) > . (B.38) 

We will find the exact definition of starting from the 5-matrix interpretation of the theory. 

It was noted previously that the turning-point field ^((Too) may be arbitrary. We will suppose that 
on the remote there are only free, on the mass-shell, particles. Formally it follows from (B.35 - 



B.37). This assumption is natural also in the S'-matrix framework [M. In other respects the choice of 



boundary condition is arbitrary. 

Therefore, we wish to describe the evolution of the system in a background field of mass-shell 
particles. The restrictions connected with energy-momentum conservation laws will be taken into 
account and in other respects background particles are free. Then our derivation is the same as in pT| . 
Here we restrict ourselves mentioning only the main quantitative points. 

Calculating the product a nm a* nm we describe a time ordered processes of particle creation and 

absorption described by D_| and D |_. In presence of the background particles, this time-ordered 

picture is slurred over because of the possibility to absorb particles before their creation appears. 
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The processes of creation and absorption are described in vacuum by the product of operators 
4>+<j)- and 4>-<j)+. We can derive (see also the generalizations of (B.21). The presence of the 
background particles will lead to the following generating functional: 



R, 



-iN( 



cp 



:), 



(B.39) 



where Ro((f>±) is the generating functional for the vacuum case, see ( B.30 ). The operator 

describes the external particles environment. 

The operator 4>*(q) can be considered as the creation and (f>i(q) as the annihilation operator and 
the product (j)*(q)(j)j(q) acts as the activity operator. So, in the expansion of N((j)*(j)j) we can leave 
only the first nontrivial term: 



Nf 



d£i(q)<f>i(<i) n ij<i>i(<i)> 



(B.40) 



since no special correlation among background particles should be expected. If the external (non- 
dynamical) correlations are present then the higher powers of 4>*4>j will appear in expansion (B.40). 
Following the interpretation of </>*</>j, we conclude that ntj is the mean multiplicity of background 
particles. 

Computing p cp we must conserve the translation invariance of amplitudes in the background field. 
Then, to take into account the energy-momentum conservation laws one should adjust to each vertex 
of in-going a nm particles the factor e _lQl9 / 2 and for each out-going particle we have correspondingly 

e ~-ia 2 q/2 

So, the product e ~ ta kq/2 e ~ia j q/2 can ^ e interpreted as the probability factor of the one-particle 
(creation+annihilation) process. The n-particles (creation + annihilation) process probability is the 
simple product of these factors if there is no special correlations among background particles. This 
interpretation is evident in the CM frame otk — (— «/?fc,0). 

After these preliminaries, it is not too hard to find that in the CM frame we have: 



n ++(qo) = n — (q ) 



e: 



90 '-' 



gr+gj ; 



e~ < 

= n(\q 



9 



(B.41) 



Computing n^ for i ^ j we must take into account the presence of one more particle: 



n+-(qo) = 0(qo) 



_ £ 1 ±h 

E n = i ne 2 <?°™ 

EOO 
n- 



gi+gi 



(fun 



&(-qo 



E w =o ne 2 q 
E c 



6(go)(l + h(q Q Pi)) + e(-q )h(-q Pi) 



and 



ri- + (q ) = Q{qo)n{q (3 2 ) + 6(-<fo)(l + n(-q p 2 )). 
Using ( |B.4l| ), JB.42|) and JB.43|) , and the definition (|B.35D we find the Green functions: 



where 



G itj {x-x',{P)) 
iGd(q, (/?)) -- 



d 4 q 



(B.42) 
(B.43) 

(B.44) 







lo xtJ 2n / 'n— 2— M) n((3 2 \q \)a + (P 2 ) 



(B.45) 
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and 



a± (/3) = _ e f CI«o|±«o)_ 



The corresponding generating functional has the standard form: 



(B.46) 



(B.47) 



where the summation over repeated indexes is assumed. 

Inserting (B.47) in the equation of state (2.S) we can hnd that (3\ — P2 — /3(E). If (3(E) is a 'good' 
parameter then Gij(x — x';(3) coincides with the Green functions of the real-time finite-temperature 
field theory and the KMS boundary condition: 



G+-(t - t') = G-+(t -t'- i/3), G-+(t - t') = G+-(t -t' + i/3), 



(B. 



is restored. The eq.(B.48) can be deduced from (B.45) by direct calculations. 



C Appendix. Local temperatures 

We start this consideration from the assumption that the temperature fluctuations are large scale. 
We can assume that the temperature is a 'good' parameter in a cell the dimension of which is much 
smaller then the fluctuation scale of temperature. (The 'good' parameter means that the corresponding 
fluctuations are Gaussian.) 

Let us divide the remote hypersurface Coo on a iV c and let us propose that we can measure the 
energy and momentum of groups of in- and out-going particles in each cell. The 4-dimension of cells 
can not be arbitrary small because of the quantum uncertainty principle. 

To describe this situation we decompose the (5-function of the initial state constraint 
product of (N c + 1) 5- functions: 

m . JV C m„ N c 

s^(p-j2ik)= / nw*-Eft.")^ (4) ^-E^). 

k=l J v=\ k=l v=\ 

where qk,u are the momentum of /c-th in-going particle in the z^-th cell and Q v is the total 4-momenta 
of n v in-going particles in this cell, v — 1, 2, N c . Therefore, 

N m u 
v=\ k=l 

The same decomposition will be used for the second (5-function of outgoing particle constraints. 
We must take into account the multinomial character of particle decomposition on TV groups. This will 
give the coefficient: 

. N , N 

— j — ' ;8 K (n - V n v ) : : -8 K (m - V m„), 

where 8k is the Kronecker's symbol. The summation over 

{ni,n 2 , ■■■,n Nc } = {n} Nc , {mi,m 2 , ...,toat c } = {m} Nc 

is assumed. 

In result, the quantity 

^v c (^Q) = £ { „, m}Nc /K™l 2 x 
x n£i {^%tH^ (4) (Q, - Er=! ik,u)^^5^(p v - TTkUvk,,)} (c.i) 
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defines the probability to find in the v-th cell the fluxes of in-going particles with total 4-momentum 
Q v and of out-going particles with the total 4-momentum P v . The sequence of these two measurements 
is not fixed. 



The Fourier transformation of (5-functions in (C.l) gives: 



2,„) 



PN c (a), 



where 

PN c (a) = PN c (ai,l,Oii,2---, <Xl,N c ; &2,1, a 2,2, «2,iv) 

has the form: 



= / n In ^r^e-— n ^^e-— A \ anm ?. ^ 

J 17, IHv , i ll V- 

v=\ \k=\ k=l ) 



Inserting 

m n 

a nm (p,q) = (-i) n+m J] k<Uv>) II 4>*iPh,u)Z{-^ 
k=l k=l 

into ( |C.2| ) we find: 



PN B {a) = exp{i£)jV l / rfa;d^'[^ + (a;)D + _(a;-a;';Q!2, I /)0_(a;')- 

-0_(x) J D-+(a; - x'- ai , v )h{x')])p^), (C.3) 

where -D-i (a; — a; ; a), and I? |_(ai — a;'; a) are the positive and negative frequency correlation functions. 

We must integrate over sets {Q}n c an d {P}n c if the distribution of momenta over cells is not fixed. 
In the result. 



R(P) = J D 4 ai (P)D 4 a 2 {P) PNc {a), (C.4) 



where the differential measure 



N c , 4 



1 = 1 (2.) 

takes into account the energy-momentum conservation laws: 

K(P,{a} Nc )= / H^Q^^^S^iP -J2Qu) 
The explicit integration gives that 

K(P,{a} Nc )~]j6W{a-a v ), 
where a is 3- vector of the CM frame. Choosing CM frame, a — (— if}, 0), 

{(3} Nc ) = / TT rf£„e£^i " E "S(E - V £?„). 

In this frame 

PWc (P)= / D(3 x {E)Dfh{E)p N M, 
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where 

D(3{E) = T[^K[E,{f}} Nc ) 



2m 

v=l 



and /3at c (/3) was denned in (C.3) with au.u = ( — i/3fc ; ^, 0) , Re/3k,u > 0, k = 1,2. 

We will calculate integrals over Pk using the stationary phase method. The equations for the most 
probable values of fit'- 

.lnK(E,{p} Nc ) = -— ln PjVc (/3), fc = l,2, (C.5) 



always have unique positive solutions /3£ V {E). We propose that the fluctuations of Pk near [3^ v are 
small, i.e. are Gaussian. This is the basis of the local-equilibrium hypothesis |)7|. In this case 1/ v 
is the temperature in the initial state in the measurement cell v and l/^u ^ s the temperature of the 
final state in the ^-th measurement cell. 



The last formulation (C.4) implies that the 4-momenta {Q}n c and {P}n c cannot be measured. 
It is possible to consider another formulation also. For instance, we can suppose that the initial set 
{Q}n c is fixed (measured) but {P}n c is not. In this case we will have a mixed experiment: (31 v is 
defined by the equation: 

E » = -dL inpN * 



and @2 v i s defined by the second equation in (C.5) 



Considering the continuum limit, N c — * oo, the dimension of the cells tends to zero. In this case 
we are forced by quantum uncertainty principle to assume that the 4-momenta sets {Q} and {P} are 
not fixed. This formulation becomes pure thermodynamical: we must assume that just {Pi} and {P2} 
are measurable quantities. For instance, we can fix {Pi} and try to find {^2} as a function of the total 
energy E and the functional of {Pi}. In this case eqs.( |C.5D become the functional equations. 

In the considered microcanonical description, the finiteness of temperature does not touch the 



quantization mechanism. Indeed, one can see from (C.3) that all thermodynamical information is 
confined in the operator exponent 



the expansion of which describes the environment, and the 'mechanical' perturbations arc described by 
the functional po(4>). This factorization was achieved by the introduction of the auxiliary field 4> and is 
independent from the choice of boundary conditions, i.e. from the choice of the systems environment. 



Wigner functions 



We will adopt the Wigner functions formalism in the Carruthers-Zachariazen formulation |39| ■ For 
the sake of generality, the m into n particles transition will be considered. This will allow the inclusion 
of the heavy ion- ion collisions. 

In the previous section, the generating functional Pn c (P) was calculated by means of dividing the 
'measuring device' on the remote hypersurface (Too into N c cells 

p Nc (a) = e- iN ^^p ( ( j>), (C.6) 

where 

N(<P;P,z) = [j:^iJdxdx' x 
x ((i> + (x)D + ^(x- x';P 2 ,v,z 2 )<t>-(x') - 

-4>-{x)D_ + {x-x';p hv ,z 1 )4>+^'))} (C7) 
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is the particle number operator. The frequency correlation functions and arc defined by 

equalities: 



D + _(x-x';P 2 ,u,z 2 ) = -il dn 1 (q)e l ^ x -^e-^^z 2 (q 2 ^) (C.8) 
D_ + (x - x', /V, zi) = i I d^e-^^-^e-^^zxiq^) (C.9) 



It was assumed that the dimension of the device cells tends to zero (N c — > oo). Now we wish to specify 
the cells coordinates. In the result we will get to the Wigner function formalism. 
Let us introduce Wigner variables |9q] : 

x — x = r, x + x — 2y : x — y + r/2, x = y — r/2. (C.10) 

Then 

N( ( j>;(3,z) = -iY:^ 1 fdn(q)drx 
x U+{y + r/2)4>-(y - r/2)z 2 (q 2 ^)e^ r e-^^ + 

+ 4>-(y + r/2)j> + {y - r /2)z 1 {q 1 ^)e- l ^ r e~^<^ dy (C.ll) 

The Boltzmann factor, e~^ i ' v ' L ^ qi - v \ can be interpreted as the probability to find a particle with the 
energy e(qi tU ) in the final (i — 2) or initial (i — 1) state. The total probability, i.e. the process of 
creation and further absorption of n particles, is defined by multiplication of these factors. Besides 
e* 92 " 1 * is the out-going particle momentum measured in the ^-th cell. 

Generally it is impossible to adjust the 4-index of cell v with coordinate y. For this reason the 
summation over v and the integration over r are performed in ( p. 11 ) independently. But let us assume 
that the 4-dimension of the cell L is higher then the scale of the characteristic quantum fluctuations 

Lq, 

L»L q . (C.12) 



One can divide the 4-dimensional y space into the L dimensional cells. Then, because of (C.12), the 
quantum fluctuations can not take away particles from this cell. Then we can adjust the index of the 
measurements cell with the index of the y space cell. 
As a result, 

N((f>; (3, z) = -ifdyf dfli(q)dr x 
x (j>+{y + r/2)0-(y - r/2)z 2 (q 2 , y y^ r e-^<^ + 

+$-(y + r/2)$ + (y-r/2)z 1 (q 1 ,y))e-^ r e-^y^\ (C.13) 

where 

/ dy = V / dy, (C.14) 

and C{v) is the dimension L of the y space cell with index v. Notice that the momentum q did not 
carry the index v (or the index y of the space cell). 

Our formalism allows the introduction of more general 'closed-path' boundary conditions. The 
presence of external black-body radiation will only reorganize the differential operator exp{N((f>*(f>j)} 
and a new generating functional p cp has the same form: 



Pcp (/3,z) = e- 4N ^ z VoW>). 
) is strictly the same as in 1 

/ drdyfo(r + y/2)n t3 (y)^{r-y/2), 



The calculation of operator N(tfi*(fij) is strictly the same as in Appendix B. Introducing the cells we 
will find that 
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where the occupation number carries the cell index y: 

n>ij( r , v) = J ^li(q)e iqr ni j (y,q) 

and (q = e(q)) 

n ++ (y,q )=n—(y 1 q )=n{y 1 (ft + f3 2 )\q \/2) = e(Pl+0a)l y )M/2 _ ^ 

n+-(y,qo) = 9(Q , o)(l + ra(y,^ 2 g , o)) + B(-g )n(y, -Ago), 
«-+(?/, <7o) = n+_(y,-g )- 

For simplicity the CM system was used. Other calculations are the same as the constant temperature 
case. 



D Appendix. Multiperipheral kinematics 

First of all two light-like 4-momenta 

Pi : 2 = P\.i - P2,\m 2 /s 

are introduced. Here are momenta of colliding particles. The final state particles momenta have 
the following representation: 

Pi = "1^2 + P'lPl + P' 1± , P2 = a 2P2 + (3' 2 Pl +P2±> 

fa = a t p2 + PiPi + kij_. (D.l) 

Sudakov's parameters, a, /?, are not independent. The mass shell conditions and the energy-momentum 
conservation laws give: 

sa'^i = to 2 + (;p' 1± ) 2 = Ef ± , sa' 2 P' 2 = E 2 1± , satfi = E? ± , 

a' 1 +a' 2 +Eai = l, ^i+^+EA = l, (D.2) 

where E$± is the transverse energy. 

We have for the multiperipheral kinematics: 

9 

Tfi 

1 « 1 » fa » ... » (3 n » & , 

S 

9 

771 

<< as << a-i << ... << a n « a 2 ~ 1 (D.3) 

s 

and the transverse momenta are restricted: 

\p' i± \ ~ \k i± \ ~ m. (D.4) 
It corresponds to small production angles in the considered CM frame: 

„ I *ii I 



, » M, (D.5) 



if the particle moves along Pi, and a similar expression exists for particles moving in the opposite 
direction, where |/3j| << \at\. In the 'central region' of the CM frame \j3i\ ~ |ct;| ~ (Ei±/E) << 1 
the angles of produced particles are large and energies are small. It should be underlined that all this 
excludes the (mini)jets formation. 
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The final-state particles phase space volume element is 

rtrr - (2" a ) 2+ " ^n d 2 gi d 2 q 2 d 2 q n + 1 1 

a °2^2+n - 167r i„ q ( +m 2 ( qi -q 2 y 2 +X 2 ' ' ' fe-^+i^+A^ q'i + 1 +\ 2 X 



x^e( a2 -a 1 )---tfnL + 1 1 (^) * =^^^n, (D.6) 
where Cy = 3 and the 4-momentum of produced particle 

h = (at - a i+ i)p 2 + (fii - fii+i)pi + [q% - qi+i)± = -a l+ ip 2 + fiiPi + [q% - qi+i)±- 
The square of pairs invariant mass: 

si = ( P [ + h) 2 = s\a 2 \, s n+1 = (k n +p' 2 f = ^k, Sl = El ai+1 



V 0-1)-L 



Q-i-l 



The energy conservation law takes the form: 

s\s 2 ■■■s n = sE 2 ± ■ --El ± . 

The trajectory of reggeized gluon is 

a(q 2 ) = 



2vr 2 J (fc 2 + A 2 )((q-fc) 2 + A 2 )' 

where A is the gluon 'mass'. If this virtuality is large, A >> m, then the gluon decay creating a pQCD 
jet. But the constraint on the multiperipheral kinematics prevent this possibility. 

Deep inelastic reactions 

For the pure deep inelastic case, when one of the initial hadrons is scattered at the angle 9 have the 
energy E' in the cms of beams whereas the another is scattered at small angle and the large transfer 
momentum Q = AEE' sin 2 (0/2) >> m 2 , is distributed to the some number of the emitted particles 
due to evolution mechanism we have |9l) (8 is small): 

x/^AefPl^)...?^), P( z) = 2 ^, (D.7) 

where 9^' = 0(0j+i — 0i) and the emission angle 9i — \ki\/(Emax(ai.0i)). 

Large angle production 

For the large-angle particle production process the differential cross section (as well as the total 
one) fall down with CM energy y/s. Let us consider for definiteness the process of annihilation of 
electron-positron pairs to n photons p9| : 

da DL = *lf dFnj 

dF n = riLi dx i dy l Q(x i - yi)@{y t - y J -i)e(y i )6(a; 4 ) x 

xQ{p-x n )Q(p-y n ), y i =]nj^, p = ln^, (D.8) 

The similar formulae takes place for subprocess of quark-antiquark annihilation into n large-angle 
moving gluons. 

At the end, one can consider the following possibilities: 

a) Pomeron regime (P); 

b) Evolution regime (DIS); 

c) Double logarithmic regime (DL); 

d) DIS+P regime; 
c) P+DL+P regime. 

The description of every regime may be performed in terms of effective ladder-type Feynman diagrams. 
This can be done using the blocks dZ n , dD n and dF n . 
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E Appendix. Reggeon diagram technique for gen- 
erating function 

We will consider, see ( |3.1l| ) 

V{q,u>;z) = <%e-"tp(q,Z;z)= , g = ln( S /m 2 )- (E.l) 

Jo w + "off + wW 

as the 'propagator of the cut Pomeron'. It will be assumed also that 

ipo(z) = -A + (1 - z)n , n > 0. 

So, the resonance short range correlations will be ignored in this definition or propagator. It was 
assumed also that the 'bare' slope a' is z independent. 



It should be underlined that the 'propagator' (E.l) is written phenomenologically. It absorbs the 



assumptions that (i) the diffraction cone shrinks with energy and (ii) the inclusive cross sections are 
universal, see (^) 

The set of principal rules concerning multiperipheral kinematics of Feynman diagrams is given in 
Appendix D. The reggeon calculus supposes that the virtuality of each line of the Feynman diagram is 
restricted. This ignores 'hard jets', later known as the pQCD jets. 

Then the v Pomeron exchange eikonal diagram has only (y + 1) ways of being cut. If the cut line 
goes through /x Pomerons, then the corresponding contributions are: 



dSl v {M^q 1 ,...,q v )) 2 y^ilV{g il UHiz) ]J V{q u ^z = \\ (E.2) 

J 1=1 i=/J,+l 

where Mj?(qi, ■ ■■,q v ) is the 'vertex function', the combinatorial coefficient is 

v „ _ (-1) ( "-" ) 2'V! 
»\{ V -H)\ 

and the phase space element is 

^=n^^-E^ 2 ( 5 -E«). (e.3) 

i=l ^ ' 1 = 1 i=l 



As usual, the contribution (|E.2|) leads to the following mean multiplicity of produced particles: 



fj,fi(s). (E.4) 



Therefore, 

ju~n(s) (E.5) 

is essential in the VHM region, where n ~ n(s) 2 is assumed. 
The impact parameter representation: 

= J <m(JLy ft^M (e.6) 



would be useful also. The contribution (E.2) describes interactions with impact parameter 

< b 2 >~ Aa'\TL(s/m 2 )/v. (E.7) 

Notice < b 2 > is the number of cut pomerons independent of /i. But, remembering that n > v and 

that the Regge model is only able to describe large distance interactions, m 2 < b 2 >> 1, one can 
conclude that the Regge pole description is valid only for 

n<n(sf. (E.8) 

Thus is why the VHM region is defined by n(s) 2 . 
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F Appendix. Pomeron with A > 



Then the cut Pomeron propagator in the impact parameter representation 

~g(b,t;z)=g(h,Oe^-^\ (F.l) 

where 

ff(b,£) = ^e A « e - b2 /^ (F.2) 

is the uncut Pomeron profile function. Using this definition one can find that the contribution of the 
cikonal diagrams gives a contribution: 

T (b, f ; z) = (l - e~ x2 ^) - 1 (l - e 2A 2 9 (b, ? )(^- 1)fi(s) -D) ( (F 3) 

where A is a constant. 

First bracket is essential for b 2 < 4a' A£ 2 . So, with exponential accuracy, the first term is equal to 

9(4a'A£ 2 - b 2 ). 
Let us now consider the second bracket. For z < 1 

b 2 < 4a' AC 2 jl + i In (l - e (*-i)«W^ J = 4a'A<e 2 7(£, z) (F.4) 

are essential. It is not hard to see that 7(£, z) decreases if z — > 1 and 7(£, z) is equal to zero for z = 1. 
In this case, by the definition of the generating function, the integral over b of J-o(h, £; z = 1) defines 
the contribution to the total cross section. So, the model predicts the production of particles in the 
ring: 

4a'£ 2 7 (£, z) < b 2 < 4a'£ 2 A. (F.5) 
if z < 1, i.e. if n < n(s). Notice also that 7(£, z = 0) = A + O^^). Then, 

^o(b,e;z = 0) = ^o(b,C;« = l). (F.6) 

So, the elastic part of the total cross section is half of the total cross section. This means, using the 
optical analogy, that the scattering on the absolutely black disk is well described. 

The last conclusion means that the interaction radii should increase with n in the VHM region. 
Indeed, as follows from ( JF. 3| ) , 

^(b.^^^i^ 2 ^)^- 1 ^'- 1 )-!), (F.7) 

at z > 1 and, therefore, 

< b 2 < B 2 = 4a'£(A£ + (z - l)n(s)) (F.8) 

are essential. 

G Appendix. Dual resonance model of VHM events 

Our purpose is to investigate the role of the exponential spectrum ( |3.19[ ) in the asymptotic region over 
multiplicity n. In this case one can validate heavy resonance creation and such formulation of the 
problem have definite advantages. 

(i) If creation of heavy resonances at n — > oo is expected, then one can neglect the dependence on 
the resonance momentum qi . So, the 'low-temperature' expansion is valid in the VHM region. 

(ii) Having the big parameter n, one can construct the perturbations expanding over l/n. 
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(iii) We will be able to show at the end the range of applicability of these assumptions. 
For this purpose, the following formal phenomena will be used. The grand partition function 



T(z,s) = ^V l cr„(s), T(l, s) = atot(s), n<y/s/m = n max (s), (G.l) 



will be introduced, see ( p.3q ). Then the inverse Mellin transformation, see (2.37) 



will be performed expanding it in the vicinity of the solution z c > of the equation of state, see (2.38): 

d 

n = z—lnT(z,s). (G.3) 

It is assumed, and this should be confirmed at the end, that the fluctuations in the vicinity of z c are 
Gaussian. 

It is natural at first glance to consider z c = z c (n, s) as an increasing function of n. Indeed, this 
immediately follows from the positivity of o~ n (s) and the finiteness of n max {s) at finite s. But one can 
consider the 'thermodynamical limit', see Sec. 2. 3.1, or the limit mo — ► 0. Theoretically, the last one is 
right because of the PCAC hypotheses and nothing should happen if the pion mass too — > 0. In this 
sense T{z, s) may be considered as the whole function of z. Then, z c = z c {n, s) would be an increasing 
function of n if and only if T(z, s) is a regular function at z = 1. 

The proof of this statement is as follows. We should conclude, as follows from eq.(G.3), that 

z c (n, s) — ► z s at n — ► oo, and at s = const, (G.4) 

i.e. the singularity point z s attracts z c in asymptotics over n. If z s = 1, then (z c — z s ) — > +0, when 
n tends to infinity |50|. The concrete realization of this possibility is shown in Sec. 3. 3. But if z s > 1, 
then (z c — z s ) — > —0 in VHM region, see Sees. 3.1, 3.5. 
On may use the estimation, see also ( 2.39 ): 

- iln-^v =lnz c (n,s)+0(l/n), (G.5) 
n cr tot (s) 



where z c is the smallest solution of ( p.3| ). It should be underlined that this estimation is independent 
of the character of singularity, i.e. the position z s is only important with 0(l/n) accuracy. 

Partition function 



Introducing the 'grand partition function' (|G.l|) the 'two-level' description means that 



xN k (q 1 ,q 2 ,...,q k ;P) = -f3F{z,s), 



(G.6) 



where e(qi) = (qf + ml) 1 / 2 . This is our group decomposition. The quantity £(q, z) may be considered 
as the local activity. So, 

d ' ~<rtotBi{n) (G.7) 



If the resonance decay forms a group of particles with total 4-momentum g, then B\ (q) is the mean 
number of such groups. The second derivative gives: 



S 2 T 



^(qi,z)S^(q 2 ,z) 



crtot{B2{qi,q2) - B 1 (qi)Bi(q 2 )} = a to tK 2 (qi,q2) 



(G.S 
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where ^2(91,92) is the two groups correlation function, and so on. One can consider B k as Mayer's 
group coefficients, see Sec. 2. 3. 3. 

The Lagrange multiplier (3 was introduced in ( |G.6| ) to each resonance: the Boltzmann exponent 
exp{— /3e} takes into account the energy conservation law J^- ej = E 1 where E is the total energy of 
colliding particles, 2E = -^i in the CM frame. This conservation law means that (3 is defined by 
equation: 

\/i= AlnT(z,/3). (G.9) 



So, to define the state one should solve two equations of state (G.3) and (G.9). 

The solution (3 C of the eq.(G.9) has the meaning of inverse temperature of the gas of resonances if 
and only if the fluctuations in vicinity of (3 C are Gaussian, see Sec. 2. 2. 2. 

On the second level, we should describe the resonance decay into hadrons. Using ( 3.24 ) we can 
write in the vicinity of z = 1: 



z — ' m 



(G.10) 



The assumptions B and D, see ( 3.21 ), were used here 
So, 

- f3F(z,s) = / n{d™>i£(mi,z)}B k (m;/3), 

k J i=l 

where m = (mi, m,2, rrik) £ was defined in (G.10) and 

B k (m;l3) = J ' J[{dSl k (q)e-^^} B k (m;q). 

i—l 

Assuming now that |qi| << m are essential, 

k 



i=l 



2rrii 



Following the duality assumption, one may write: 



B k (m) = 4(m)nKe Am '} 



and B k (m) is a slowly varying function of m = = (mi, m.2, m k ): 

N k (m) ~ 6 fc 

In the result, the low-temperature expansion looks as follows: 

2 k l 2 m%{ 9 R ) k b k 



(3T{z, s) = J2 



f]3k/2 



^ mm 7+3/2 e ( 2 -l)« R (m)-(/3-/3 )m 



(G.ll) 



(G.12) 



(G.13) 



(G.14) 



(G.15) 



We should assume that (/3 — (3q) > 0. In this sense one may consider l//3o as the limiting temperature 
and the above mentioned constraint means that the resonance energies should be high enough. 



Thermodynamical parameters 



50 



Remembering that the position of the singularity over z is essential, let us assume that the resonance 
interactions can not renormalize it, i.e. that the sum ( G.15| ) is convergent. Then, leaving the first term 
in the sum ( G.15| ), 

R/~i poo 

-(3f(z, S ) = m °°, 1 / dm(m/m r+ 3 / 2 e^- 1 ^ n ^-^-M m . (G.16) 

P i/2 J, n 

We expect that this assumption is satisfied if 

Jm bifil 6 '-* 

for 

nm _ n , 
n — > oo, s — > oo, — -j=r = << 1 (G.18) 

\J S Tlmax 

So, we would solve our equations of state with the following free energy: 

= ^ r d (—)(—y'^e- A ^\ (G.19) 



/3 3/2 Jm m o rn 



where, using ( |3.20D , 

7' = 7 + 2(z - 1)uq + 5/2 = 2(z - l)n* A = m (/3 - A>) > Q, a = const. (G.20) 
We have in terms of these new variables the following equation for z, 

2an$ d r( 7 ',A) 



03/2 5 y A7 ' • 

The equation for (5 takes the form: 

_ amp r(y + 1, A) 

fimra ~~ ^3/2 At' +1 
where n max = (^/s/mo) and r(A, 7') is the incomplete T-function: 

dxx 1 ~ x e~ x . 



(G.21) 



(G.22) 



roc 

r(7',A)= / 



Asymptotic solutions 



Following physical intuition, one should expect the cooling of the system when n — > 00, for fixed 
y^s, and heating when n max — ► 00, for fixed n. But, as was mentioned above, since the solution of 
eq.( G.22 ) (5 C is defined by the value of the total energy, one should expect that (5 C decreases in both 
cases. So, the solution 

<9A <9A 

A c > 0, — — < at n -> 00, —2 < at s -> 00 (G.23) 
on os 

is natural for our consideration. 

The physical meaning of z is activity. It defines at (5 = const the work needed for one particle 
creation. Then, if the system is stable and T(z, s) may be singular at z > 1 only, 

dz d z 

— - > at n -> co, — - < at s -> 00. (G.24) 
an os 
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One should assume solving equations ( G.2l[ ) and ( |G.22 ) that 

z^ 1 ^ 1 '^ «r(y c +l,A c ). (G.25) 
die A? c 

This condition contains the physical requirement that n << n max . In the opposite case, the finiteness 
of the phase space for mo ^ should be taken into account. 

As was mentioned above, the singularity z s attracts z c at n — > oo. For this reason one may consider 
the following solutions. 

A. z s = oo: z c » A, A << 1. 

In this case 

A- 7 'r( 7 ',A) ~e 7 ' ln(7 ' /A) . (G.26) 
This estimation gives the following equations: 

n = Ci7'ln( 7 7A)e 7 ' ln(7 ' /A) ,^— = C 2 A 7 'ln(^) « 1, (G.27) 



where C{ — 0(1) are the unimportant constants. The inequality is a consequence of flaH). 
These equations have the following solutions: 

<< 1, i c ~ Inn >> 1. (G.28) 



In? 



Using them one can see from ( G.5 ) that it gives 

a n < 0(e~ n ). (G.29) 

B. z s = +1: z c -► 1, A c << 1. 

One should estimate r(Y, A) near the singularity at z = 1 and in the vicinity of A = to consider 
the consequence of this solution. Expanding T(j' , A) over A at 7' — > 0, 

r( 7 ', A) = T( 7 ') - A~>'e- A + 0(At' +1 ) ~ — + 0(1). (G.30) 

7 

This gives the following equations for 7 ': 

n = c , Yln(l/A)-l eVln(1/A) _ 

7 

The equation for A has the form: 

n ma:c = C^' +1 ) ln ( 1 / A >. (G.32) 

Where C[ = 0{1) are unimportant constants. 
At 

< 7 'ln(l/A) - 1 « 1, i.e. at ln(l/A) « n « ln 2 (l/A), (G.33) 

we find: 

^Riky- (G - 34) 

Inserting this solution into ( p.32 ): 

A c — . (G.35) 

Umax 

It is remarkable that A c in the leading approximation is n independent. By this reason 7 £ becomes n 
independent also: 

This means that 

<r„ = 0(e~ n ) (G.37) 
and obeys the KNO scaling with mean multiplicity n = Hq ]n(n max ) . 
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H Appendix. Correlation functions in DIS kine- 
matics 

Considering particle creation in the DIS processes, one should distinguish correlation of particles in 
the (mini)jets and the correlations between (mini)jets. We will start from the description of the jet 
correlations. One should introduce the inclusive cross section for the v jets creation 

<S> ( J )v (k 1 ,k 2 ,.... 7 k u ;q 2 ,x), 

where fcj, i = 1,2, ...,n are the jets 4-momentum in the DIS kinetics, — q 2 >> A 2 . Having <& v we can 
find the correlation functions 

N^ r \ki 1 k 2 ,....,k 1J ]q 2 1 x) 1 

where (r) = n, ...,r u and r.i = {q,g) defines the sort of created colour particle. It is useful to introduce 
the generating functional 

~ n 

F ab (q 2 ,x;w) = Y, / dSl n {k) J[ w r ' (h) \a a n b (k u k 2 , k n ; q 2 , x)\ 2 , (H.l) 

n J i=l 

where a° 6 is the amplitude, dfl n (k) is the phase space volume and w Ti {ki) are the arbitrary functions. 
It is evident that 

F ab (q 2 ,x;w) \ w=1 = D ab (q 2 ,x). (H.2) 

The inclusive cross sections 

k 2 , k v ; q 2 ,x) = I] *_jr«V, x; w) \ w=1 . (H.3) 

The correlation function 

N^(k;q 2 ,x) = f[-^L-^nF ab (q 2 ,x-,w) \ w=1 . (H.4) 

We can find the partial structure functions D ab (q 2 , x; n), where n is the number of produced (time- like) 
gluons, using there definitions. 

It will be useful to introduce the Laplace transform over the variable ln(l/ic): 

F a \ q 2 ,x;w)= [ ±(iY r \q 2 , J]W ) (H.5) 
JRej<o 2vri \xj 

The expansion parameter of our problem a s ln(— q 2 /X 2 ) ~ 1. For this reason one should take into 
account all possible cuts of the ladder diagrams. So, calculating D ab (q 2 , x) in the LLA all possible cuts 
of the skeleton ladder diagrams are defined by the factor |Q: 

-{TfG r Tf\, (H.6) 

i.e. the cut line may not only get through the exact Green function G r (kf) but through the exact 
vertex functions r" b (<7i, qi+i, ki) also (qf,qf +1 are negative). We have in the LLA (see Appendix D) 

A << -tfc << -q l+1 « -q 

and 

X < Xi+i < Xi < 1. 

Following our approximation, see previous section, we would not distinguish the way in which the cut 
line goes through the Born amplitude 

af = {(rf) 2 G r }. 
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We will simply associate wTma" 6 to each rung of the ladder. 

Considering the asymptotics over n, the time-like partons virtuality ki ~ —q 2 /y% should be maxi- 
mal. Here y% is the fraction of the longitudinal momentum of the jet. Then, slightly limiting the jets 
phase space, 



Inkf = In \q l+1 \ z (1 + O(ln(l/a0/l<fc+ir))- 



(H.7) 



As a result, introducing r, = ln(qf /A 2 ), where a s (q 2 ) = 1//3t, = (UN/3) - (2n//3) in the LLA 
variable, we can find the following set of equations: 



d 

T Q^fab(q 2 ,j;w) =^2ip r ac {j)w r {T)f ab {q 2 ,x;w), 

rl dx 



where 



<PacU) = Vac(j) 



-X 3 P ac (x) 



(H. 



(H.9) 



and P ac (x) is the regular kernel of the Bethe-Salpeter equation |fT8| . At w = 1 this equation is the 
ordinary one for D ab (q 2 ,x). 

We will search the correlation functions from eq.( H.8 ) in terms of the Laplace transform 

Kb " ( fc i > fc 2 K ; q 2 , 3 ) = n ab " ( k ; <? 2 > i ) 

Let us write: 



fab(q 2 ,j;w) = d ab (q 2 , j) exp 



^ 1/ Z — 1 V 



))n%"(k;q 2 ,j) 



(RAO) 



Inserting ( H.10 ) into (H.8) and expanding over (w — 1) we find the sequence of coupled equations. 
Omitting the cumbersome calculations, we write in the LLA that 



= d aCl (3,ri)<p r c \ C2 (j)d C2C3 (j,T2) ■ ■■fZ Cu+1 (j)d Cl , +1 b(j,Tu+i) 
One should take into account the conservation laws: 

Tl ■ T 2 ■ ■ ■ T u+ l = T, Tl < T 2 < . . . < Tu+l < T. 

Computing the Laplace transform of this expression we find 

( r l.7"2, ■■■■,T„;q 2 ,x). 



(H.ll) 
(H.12) 



The kernel d a b(j,r) was introduced in (H.ll). Let us write it in the form: 



dab(j,r) =^2<? 



dab(j) t u € 

v + — V- 



where 

d qq = V s — (fgg, d qg = V s — (fqq, d qg — (fgq, d gq = Ifqg 

and 

v °=2 {vii + t Pgg+ (T [(Vqq - Vggf - ^fVggVgq] V2 } 
If x << 1, then (j — 1) << 1 are essential. In this case fT8| , 

Vgg ~ Vgq » fqg ~ <Pqq = 0(1). 

This means the gluon jets dominance and 

n g gg = <P BB + 0{l). 



(H.13) 

(H.14) 
(H.15) 

(H.16) 

(H.17) 
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One can find the following estimation of the two-jet correlation function: 



<P(n,T 2;j ,T) = O (max{(ri/r)*'", fo/r)*", (n/r 2 )^}} . (H.18) 

This correlation function is small since in the LLA t\ < ti < t. This means that the jet correlation 
becomes high if and only if the mass of the correlated jets are comparable. But this condition shrinks 
the range of integration over t and for this reason one may neglect the 'short-range' correlations among 
jets. Therefore, as follows from (H.10), 

f ab (q 2 ,j;w) = d gg {Tj)exp^cp gg J ^- w f(r')| (H.19) 

We will use this expression to find the multiplicity distribution in the DIS domain. 

Generating function 



To describe particle production, one should replace: 

w r lmaf -> w r Jmaf , 
where is the probability of n particle production, 

n 

Having v jets, one should take into account the conservation condition n\ + ni + 
reason, the generating functions formalism is useful. In the result, one can find that if we take (H.19) 

w 9 = w 9 (t, z), w 9 {t, z) \ z=1 = 1, (H.21) 



(H.20) 
n v = n. For this 



then fab{q 2 j',w) defined by (H.19) is the generating functional of the multiplicity distribution in the 
l j representation'. In this expression w 9 (t, z) is the generating function of the multiplicity distribution 
in the jet of mass |fc| = Ae 7 "/ 2 . 
In the result, see (H.5), 



F ab {q 2 ,x-w) cx 



Rej< 



2iri 



[1/xYe 



where 



rT dr' 

u(t,z) = I —w 9 (t',z). 



(H.22) 



(H.23) 



Noting the normalization condition (H.21), 

uj(t,z=1) = ]iit. (H.24) 

The integral flH.22j ) may be calculated by the steepest descent method. It is not hard to see that 

j ~ j c = 1 + {ANlu{t, z)j ln(l/o;)} 1/2 (H.25) 



is essential. Notice that j — 1 << 1 should be essential but we find, instead of the constraint ( p.26D , 
that 

uj(t,z) « lu{l/x). (H.26) 



In the frame of this constraint, 



F 



a (q , x; w) cx exp z) ln(l/x)| 



(H.27) 
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Generally speaking, there exist such values of z that j c — 1 ~ 1. This is possible if lu(t,z) is a 
regular function of z at z — 1. Then z c should be an increasing function of n and consequently oj(t, z c ) 
would be an increasing function of n. Therefore, one may expect that in the VHM domain j c — 1 ~ 1. 

Then j ~ 1 + u;(t, z)/ ln(l/a;) would be essential in the integral (H.22). This leads to the following 
estimation: 

F ab (q\x;w) oc e~^ T ' z l 

But this is impossible since F ab (q 2 ,x;w) should be an increasing function of z. This shows that the 
estimation (H.27) has a finite range of validity. 

Solution of this problem with unitarity is evident. One should take into account correlations among 



jets considering the expansion ( H.1C ). Indeed, smallness of n 
of Yl" w Ti (Ti,z) in the VHM domain. 



may be compensated by large values 



I Appendix. Solution of the jets evolution equation 

One may neglect quark jets in the VHM region since the gluons mean multiplicity n g > n q - quarks 
multiplicity pl| , |S5| and in the VHM region the leftmost singularities are important. Then we can 
write @: 

J^-Ct, z) = ^T,(r, z) £ dr\T 3 {r\ z) - 1), (LI) 

where r = ln(q 2 /A 2 ) and Tj(r,z) is the generating function of the distribution over the number of 
gluons w n (tau): 

T 3 (t,z) = ^z n W „M, T j (r,z= 1) - 1. (1.2) 

n 

We search a solution in the VHM region, where 

n >> rij oc expjv^ar}, a = 12/11. (1.3) 
Let us consider the following solution: 



It is useful to introduce 



w n = [ ^Ye- an /^. (1.4) 



a k (r) = -Y,^ 1 ^(r) (1.5) 

' 71=1 



for this solution. Inserting (L4) in this expression, 

a k (r) =n k {T)fi kl (1.6) 

where (3% (i) should be positive and (ii) r independent. 

These conditions are satisfied for the following values of k. Indeed, at k » 1, 

The generating function Tj has the following form in terms of a^: 

OO 

T j (r,z) = J2( l ^) k Mr)- (1-8) 



Inserting ( |I.8| ) into ( |I.l[ ) and assuming that ft is a t independent quantity, we find the following 
recurrent equation for ft: 

fa = \ E Y M ~ ki ~ 2 nr Pk - (L9) 
fei=i 1 
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Therefore, if 

(1.10) 

then we can neglect last term in the right hand side of dO| ) and in this case Pk are positive and r 
independent. Noting that in ( |l~7| ) n ~ kflj are essential the inequality (LK) means that the solution 
(|hj) is correct if 

n>>fijhinj, (1. 11) 



i.e., only for this value of n w n has the form (1.4) and the corresponding generating functional has 
singularity at 

z„ = 1 + — . (1.12) 

Tin 



J Appendix. Condensation and type of asymptotics 
over multiplicity 

It is important for the VHM experiment to have an upper restriction on the asymptotics. We wish to 
show that a n decreases faster than any power of 1/n: 



a n < 0(l/n). 



(3.1) 



To prove this estimation, one should know the type of singularity at z = 1. 

One can imagine that the points, where the external particles are created, form the system. Here 
we assume that this system is in equilibrium, i.e. there is not in this system macroscopical flows of 
energy, particles, charges and so on. 

The lattice gas approximation is used to describe such a system. This description is quite general 
and does not depend on details. Motion of the gas particles leads to the necessity to sum over all 
distributions of the particles on cells. For simplicity we will assume that only one particle can occupy 
the cell. 

So, we will introduce the occupation number ct, = ±1 in the i-th cell: Oi = +1 means that we have 
no particle in the cell and cr^ = — 1 means that a particle exists in a cell. Assuming that the system 
is in equilibrium, we may use the ergodic hypothesis and sum over all 'spin' configurations of <Xj, with 
the restriction: af = 1. It is evident that this restriction introduces the interactions [100]. 

The corresponding partition function in temperature representation p2[ 



p(P,H) 



(3.2) 



where integration is performed over |<t(x)| < oo and, considering the continuum limit, Da — Y[ x da(x). 
The action 

(J.3) 
(J.4) 



S x (a) = J dx |^(V<t) 2 - uja 2 + ga A - Xa 



where 



i , Per \ Per / Per 



1/2 



PH. 



and 1/Pcr is the critical temperature. 



Unstable vacuum 



We start this consideration from the case w > 0, i.e. assuming that P > P, 
state is degenerate if H 



In this case the ground 



0. The extra term 



aH in (1.3) can be interpreted as the interaction with 

-1 and is related 



external magnetic field H . This term regulates the number of 'down' spins with a 
to the activity: 



,1/2 



J3H 



(J.5) 
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i.e. H coincides with the chemical potential. 
The potential 

v(a) = -Loa 2 + ga 4 , lo > 0, (J.6) 

has two minima at 



cr± = ±y / uj/2g. 

If the dimension d > 1, no tunnelling phenomena exist. But choosing H < the system in the correct 
minimum (it corresponds to the state without particles) becomes unstable. The system tunneling into 
the state with an absolute minimum of energy. 

The partition function p(/3, z) becomes singular at H — because of this instability. The square 
root branch point gives 

Imp(M = ^e--(«/" 2 , a t >0. (J.7) 
Note, lmp(b, z) = at H = 0. Deforming the contour in the Mellin integral over z on the branch line, 

?4 



In this integral 



z c cx exp I > (J-9) 



is essential. This leads to the following estimation: 

p„ c< e~ 3 ^ 1/3 " 2/3 <0(l/n). (J.10) 

It is useful to note at the end of this section that 

(i) The value of p n is defined by Imp(6, z) and the metastable states, the decay of which gives a 
contribution into Kep(b, z), are not important. 

(ii) It follows from (fL9|) that in the VHM domain 

H - H c - lnz c - (l/n) 1/3 -> 0. (J. 11) 

So, the calculations are performed for the 'weak' external field case, when the degeneracy is weakly 
broken. It is evident that the life time of the unstable (without particles) state is large in this ca se an d 
for this reason the semiclassical approximation is correct. This is an important consequence of ( 3.46 ). 

Stable vacuum 



Let us consider now u> < 0, i.e. (3 < (3 cr . The potential (J.6) has only one minimum at a = in 
this case. The inclusion of an external field shifts the minimum to the point a c = <r c {H). In this case 
the expansion in the vicinity of a c should be useful. As a result, 



p(fi, z) = exp{ / dx\a c - W(a c )}, (J.12) 
where W(a c ) can be expanded over a c : 

W(a e ) = Y,j [l[{dx k a B (x k ;H)}b l (x 1 ,...,x l ). (J.13) 

l J k 

In this expression, bi(x%, xi) is the one-particle irreducible Green function, i.e. bi is the virial 
coefficient. Then a c can be considered as the effective activity of the correlated /-particle group. 
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The sum in (J. 13) should be convergent and, therefore, \s c 



if \H\ 



But in this case 



the virial decomposition is equivalent to the expansion over the inverse density of particles 25 . In the 
VHM region it is high and the mean field approximation becomes correct. In result is 



|A| 
4.9 



1/3 



co if |A| — > oo, 



and 



3 A 



4/3 



2/3' 



-1/2 



We can use this expression to calculate p n . In this case 

z c cx e 4gn — > oo at n — > oo, 

is essential and in the VHM domain 

Pn oc e" 4 ^ < 0(e~ n ). 



(J.14) 



(J.15) 



(J.16) 



(J.17) 



This result is an evident consequence of vacuum stability. It should be noted once more that the 
conditions ( 3.46| ) considerably simplify calculations. 



K Appendix. New multiple production formalism 
and integrable systems 



S'-matrix unitarity constraints 



To explain our idea, let us consider the spectral representation of the one-particle amplitude: 



A 1 (x 1 ,x 2 :E) = — — , e -> +0, 

hj — hi n — IE 



(K.1) 



It describes the transition of a particle with energy E from point x\ to X2 ■ According to our general 
idea, see introduction to Sec. 2.1, we will calculate 



R\(E) = J dxidx 2 Ai(xi,x2;E)A* 1 (xi,x2;E). 



(K.2) 



The integration over the end points x\ and x 2 is performed only for the sake of simplicity. 

Inserting (K.l) into (|K.2|) and using ortho-normalizability of the wave functions ^f n (x) we find that 



eR 1 {E)=eY, n 



E-E n -ie 



-LV / 1 1 \ 

2i L-m \ E-E n -ie E-E n +is J 



= ImE„ e^W-U = * En HE - E n ). 



(K.3) 



On other hand, the closed-path amplitude, offered for calculation in [101 1 

Cx (E) = E„ / dx^^ = E„ T^-r e = 



= E„ {Pe^ + i*5(E - E n )) = E„ + ieR^E). 

So, we wish to calculate only the imaginary part of the closed-path contribution: 

eR(E) = lmd(E). 



(K.4) 
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Notice the extra factor e in the left hand side. 

The reason for this choice is evident: the real part of Ci(E) is equal to zero at E = E n , i.e. did 
not contribute to the measurable. To calculate the bound states energy spectrum, it is enough to know 
the only imaginary part of the closed-path amplitude. 

This property is not accidental. It is known as the optical theorem and is the consequence of 
the total probability conservation principles. Formal realization of it is the unitarity condition for the 
S-matrix: SS + = I. In terms of the amplitudes A, S = I + iA, the unitarity condition presents an 
infinite set of nonlinear operator equalities: 

iAA* = A - A*. (K.5) 

Notice that expressing the amplitude by the path integral one can see that the left hand side of this 
equality offers the double integral and, at the same time, the right hand side is the linear combination 
of integrals. Thus, the continuum contributions into the amplitudes should be canceled to provide the 
conservation of total probability. In this sense it is a necessary condition. 

Indeed, to see the integral form of our approach, let us use the proper-time representation: 



(K.6) 



A 1 (x 1 ,x 2 ;E) = J2^n(xi)K(^)i / dTe^ E ~ E ^ T 
and insert it into ( |K.2| ): 

/•OO 

R 1 (E)=J2 dT + dT_e-^ T + +T -^e^ E - E ^ T +~ T -\ (K.7) 

We will introduce new time variables instead of T± : 

T ± =T±t, (K.8) 

where, as follows from the Jacobian of transformation, |t| < T, < T < oo. But we can put |r| < oo 
since T ~ 1/e-^oo is essential in the integral over T. As a result, 

Pl (E) = 2nY / dTe~ 2iT (K9) 

n J J — oo 

In the last integral, the continuum of contributions with E ^ E n are canceled. Note that the product 



of amplitudes AA* was 'linearized' after the introduction of 'virtual' time [ 103 1 r = (T+ — T_)/2. 

We wish to calculate the density matrix p((3, z) including the consequence of the unitarity condition 
cancelation of unnecessary contributions. Here we demonstrate the result and the intermediate steps we 
will formulate, without proof, as the statements offered in jl^, where the formalities are described. 

Dirac measure 



The statement, see Q and references cited therein, 

SI. The unitarity condition unambiguously determines contributions in the path integrals for p 

looks like a tautology since e lS ^ x \ where S(x) is the action, is the unitary operator which shifts a 
system along the trajectory^. So, it seems evident that the unitarity condition is already included in 
the path integrals. 



The rule as the path integrals should be calculated is weel known, see e.g. [102 . Nevertheless 
the general path-integral solution contains unnecessary degrees of freedom (unobservable states with 
E ^ E n in the above example). We would define the path integrals in such a way that the condition 



n It is well known that this unitary transformation is the analogy of the tangent transformations of 
classical mechanics [103]. 
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of absence of unnecessary contributions in the final (measurable) result be loaded from the very begin- 
ning. Just in this sense, the unitarity looks like the necessary and sufficient condition unambiguously 
determining the complete set of contributions. 

S2. The to- into n-particles transition (unnormalized) probability R nm would have on the Dirac 
measure the following symmetrical form: 



R nm (pi, -,Pn,qi, -,q m ) =< IlfcLi \ T (lk\u)\ 2 Uk=i \ T (Pk;u)\ 2 > u = 

= e^M/M^w-^M nr=i |r(%;w)| 2 nLi l r fe")l 2 
^0(«)nr=iirfe; W )| 2 nLiir(p fc ;«)l 2 . 



(K.10) 



Here p(q) are the in(out)-going particle momenta. It should be underlined that this representation is 
strict and is valid for arbitrary Lagrange theory of arbitrary dimensions. The cikonal approximation 
for inelastic amplitudes was considered in [ 1 04 1 

The operator O contains three element. The Dirac measure DM, the functional U(x, e) and the 
operator K(j, e). 

The expansion over the operator 



K(j,e) = ±Re 



dxdt 



5j(x, t) Se(x, t) 



dxdtj(x, t)e{x, t) 



(K.ll) 



generates the perturbation series. We will assume that this series exist (at least in Borel sense). 
The functional U(u,e) and So(u) are defined by the equalities: 



So(u) = (S (u + e)- S (u-e))+ 2Re / dxdte{x, t)(d 2 + m 2 )u(x, t), 



U(u,e) = V(u + e)-V(u-e) 



2Re 



dxdte(x, t)v'(u), 



(K.12) 



(K.13) 



c 4 



where So (it) is the free part of the Lagrangian and V(u) describes inte ractions. The quantity So{u) is 
not equal to zero if u have nontrivial topological charge (see also |l05| ). 

According to SI, considering motion in the phase space (u,p) the measure DM(u,p) has the Dirac 
form: 



DM(u, p) = Yl du ( x > t)dp(x, t)5 (u- 



SH 3 (u,p) 
dp 



6[p + 



SHj(u,p) 
Su 



with the total Hamiltonian 



Hj{u,p)= J dx{^p 2 + ^{X7u) 2 +v(u)~ju}. 



(K.14) 



(K.15) 



This last one includes the energy ju of quantum fluctuations. 
The measure (K.14) contains following information: 
a. Only strict solutions of equations 

ii _5Hjppl =0 . + SH^_ =Q 
op ou 



(K.16) 



with j = should be taken into account. This 'rigidness' of the formalism means the absence of 
pseudo-solutions (similar to multi-instanton, or multi-kink) contribution. 

b. Pnm is described by the sum of all solutions of eq.(K.16), independently from their 'nearness' 
in the functional space; 

c. p nm did not contain the interference terms from various topologically nonequivalent contribu- 
tions. This displays the orthogonality of corresponding Hilbert spaces; 

d. The measure (K.14) includes j(x) as the external adiabatic source. Its fluctuation disturbs the 
solutions of eq.(K.16) and vice versa since the measure (K.14) is strict; 



61 



e. In the frame of the adiabaticity condition, the field disturbed by j(x) belongs to the same 
manifold (topology class) as the classical field defined by (K.16) [ 105 1 . 

f. The Dirac measure is derived for real — time processes only, i.e. (K.14) is not valid for tunneling 
ones. For this reason, the above conclusions should be taken carefully. 



g. It can be shown that theory on the measure (K.14) restores ordinary (canonical) perturbation 
theory. 

The parameter T(q; u) is connected directly with external particle energy, momentum, spin, po- 
larization, charge, etc., and is sensitive to the symmetry properties of the interacting fields system q 
For the sake of simplicity, u(x) is the real scalar field. The generalization would be evident. 

As a consequence of (K.14), T(q;u) is the function of the external particle momentum q and is a 
linear functional of u{x): 



T(q;u) 



dxe iqx 



SS (u) 
Su(x) 



dxe lqx (d 2 + m 2 ) u (x) , q 2 = m 2 



(K.17) 



for the mass m field. This parameter presents the momentum distribution of the interacting field u(x) 
on the remote hypersurface if u(x) is the regular function. Notice, the operator (d 2 + m 2 ) cancels 
the mass-shell states of u(x). 

The construction ( K.17|) means, because of the Klein-Gordon operator and the external states being 
mass-shell by definition ]40[, the solution p nrn = is possible for a particular topology (compactness 
and analytic properties) of quantum field u{x). So, T(q;u) carries the following remarkable properties: 

- it directly defines the observables, 

- it is defined by the topology of u(x), 

- it is the linear functional of the actions symmetry group element u{x). 



Notice, the space-time topology of u(x, t) becomes important in calculating integral (3.2) by parts. 
This procedure is available if and only if u(x, t) is the regular function. But the quantum fields 
are always singular. Therefore, the solution T(q;u) = is valid if and only if the quasiclassical 
approximation is exact. Just this situation is realized in the soliton sector of sin-Gordon model. 

Despite evident ambiguity T(q; u) carries the definite properties of the order parameter since the 
opposite solution p nm = can only be the dynamical display of an unbroken symmetry!* 3 ], i.e. of the 
nontrivial topology of interacting fields, as the consequence of unbroken symmetry. 

If (K.16) have nontrivial solution u c (x,t), then this 'extended objects' quantization problem [ |107| 
arises. We solve it by introducing convenient dynamical variables fl5| . The main formal difficulty, see 
e.g. [|108| , of this program consists of transformation of the path-integral measure which was solved in 

i05)pr ' 



Then 



S3. The measure (K.14) admits the transformation: 



u. 



(K.18) 



12 The following trivial analogy with ferromagnetic may be useful. So, the external magnetic field 
~ p,, if /2 is the magnetics order parameter, and the phase transition means that \x ^ 0. r(g, u) has 



H 

just the same meaning as TL. 

13 The ^-matrix was introduced 'phenomenologically', see also the example considered in [106, |95| 



postulating the LSZ reduction formulae, sec cq.(B.l). So, the formal constraints, e.g. the Haag theorem, 
would not be taken into account on the chosen level of accuracy. 



14 Number of problems of quantum mechanics was solved using also the 'time sliced' method [ 109 1 . 
This approach presents the path integral as the finite product of well defined ordinary integrals and, 
therefore, allows perform arbitrary space and space-time transformations. But transformed 'effective' 
Lagrangian gains additional term ~ h 2 . Last one crucially depends from the way as the 'slicing' 
was performed. This phenomena considerably complicates calculations and the general solution of 
this problem is unknown for us. It is evident that this method is especially effective if the quantum 
corrections ~ h play no role. Such models are well known. For instance, the Coulomb model in 
quantum mechanics, the sine-Gordon model in field theory, where the bound-state energies are exactly 
quasiclassical. 
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and the transformed measure has the form: 



DM(u, P ) = n dt(t)drj(t)S [i - 6 (rj + ) , (K.19) 

where hj(^,rj) = Hj(u c ,p c ) is the transformed Hamiltonian. 



It is evident that (£, 77) are parameters of integration of eqs.(K.16) and they form the factor space 
W = G/G c . For instance, if one particle dynamics is considered, then one may choose £ = x(0) and 
r/ = p(0). One may consider also the following possibility: 



and 



du 



y/2( V -v(u)) 



r\ = p 2 /2 + v(x) 

In this terms hj = rj — j(t)u c (£,, rf) and new Hamilton equations have the form: 

i = l-j 9 -^,r ] = j d -^. (K.20) 

So, we have at j = 0: £ = t + to and ?/ = rjQ. By this reason 
S 1 ^. The (action, angle)-type variables are mostly useful. 

According to ( |K.18| ) there exists transformation of the perturbation generating operator: 

55. The operator K has following transformed form: 

2K = y dt{j r e < +] ri -e ri } 1 (K.21) 

in the factor space, where jx,£x, X = £,17, are new auxiliary variables. 

As a result of mapping of the perturbation generating operator K on the manifold W the equations 
of motion became linearized: 

DM = J! S (i - _ 5 { fj - jv ). (K.22) 

Then 

56. If the Feynman's ie-prescription is adopted, then the Green function of eq.(K.2i) 

g{t - t') = 6(t - t') (K.23) 

Later on we will consider the soliton sector of sin-Gordon model. In this case £i is the coordinate 
and rji is the momentum of i-th soliton and N is the number of solitons. 

Expansion of exp{K(je)} gives the 'strong coupling' perturbation series. Its analysis shows that 

57. Action of the integro- differential operator O leads to the following representation: 

Rnm(P, q) = J w |d£(0) • g|oyi?« m (>, q) + 

drKO) • ^jyi^M} . (K.24) 

This means that the contributions into R n m(jp, q) are accumulated strictly on the boundary 'bifurcation 
manifold' dW |11C[ |, i.e. depends directly on the topology of W. 

Multiple production in sin-Gordon model 
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Let us consider now the completely integrable sin-Gordon model. For the sake of simplicity the 
integral: 

R 2 (q) = e -^U-,e) f DM{u lP )\T{q-u)\ 2 e lSo{u) - lU ^' e \ (K.25) 



where T(q;u) was defined in (p. 2;), will be calculated. 
The effective potential of the sin-Gordon model 

2m 2 f 

U(u]y;e c ) = — / dxdt sin Aujv (sinAe — Ae) (K.26) 

with 

du c du c 

e ' = e r-^-en- W CK.27) 

Performing the shifts in ( K.22 ) : 

m - m + f dt' 9 (t - t'^ut') = &(t) + m, 



Vi{t) -» »*(*) + / - Oi*7,i(*') = + V'M, (K.28) 



we can get the Green function g(i — i') into the operator exponent: 



K(ej) = ~ / dtdt'eCi - ' ee(i') + #(*) • ^(*0>- ( K - 2 9) 



Note the Lorentz noncovariantness of our perturbation theory with Green function (K.23). 
As a result: 

W 8h 
D N M(i,r 1 )=^^dUt)d Vl (t)S(i l -u J ( V + r ] / ))S(r ]l ), ufa) = ^ (K.30) 

with 

u w = u N (x;£ + £,\v + v')- (K.31) 

Using the definition: 

DxS(x) — J dx(0) — J dxo, 

the functional integrals on the measure ( K.30| ) are reduced to the ordinary integrals over initial data 
(£,rf)o. These integrals define zero modes volume. Notice that the zero-modes measure was defined 
without the Faddeev-Popov anzats. 

We would divide the calculations into two parts. First of all, we would consider the quasiclassical 
approximation and then we will show that this approximation is exact. 

This strategy is necessary since it seems to be important to show the role of quantum corrections 
noting that for all physically acceptable field theories R n m = in the quasiclassical approximation. 

The iV-soliton solution un depends upon 2N parameters. Half of them, N, can be considered as 
the position of solitons and the other TV as the solitons momentum. Generally at \t\ — > oo the un 
solution decomposed on the single solitons u s and on the double soliton bound states Ub [ 1 1 1 1 : 

ni n 2 

u N (x,t) = ^2u sJ (x,t) +^2u byk (x,t) +0{e~ w ) 

3 = 1 fc=l 

Note that this asymptotic is achieved if £j — ► oo or/and r)i — > oo. This last one defines the bifurcation 
line of our model. So, the one soliton u s and two-soliton bound state Ub would be the main elements 
of our formalism. Its (£, rf) parametrizations have the form: 

4 A 2 

u s (x;^,r]) = — — arctan{exp(ma;cosh/?ry — £)}, f3 — — (K.32) 
A 8 
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and 



u b (x;£,r]) 



4 . /3m mx sinh sm, cos ^ - £ 2 , 

— arctanltan 4 ■== }. 

A 2 ma;C osh^sm^a-a 



Performing last integration we find: 



A' 



A 



UN 



u N = u N (r) + r/,£o + w(t) + £')• 



uj(t) = / dt'9(t-t')u)(n + T]'){t') 



where 
and 



In the quasiclassical approximation £' = rj = we have: 

un = u N (x;i]o,£ + Lu{rj )t). 

Notice that the surface term 



Then 



dx"d^e tqx u N ) = 0. 



d 2 xe iqx (d 2 + m 2 )u N (x, t) = -(q 2 - m 2 ) / d 2 xe iqx u N (x, t) = 



(K.33) 

(K.34) 

(K.35) 
(K.36) 

(K.37) 
(K.38) 
(K.39) 



since q 2 belongs to the mass shell by definition. The condition (K.38) is satisfied for all q M ^ since 
Un belongs to Schwarz space (the periodic boundary condition for u(x,t) do not alter this conclusion). 
Therefore, in the quasiclassical approximation i? 2 = 0. 

Expanding the operator exponent in ( K.34 ) we find that action of operators ff creates terms 



d 2 xe iqx 6{t - t'){d 2 + m 2 )u N {x, t) ^ 0. 



(K.40) 



So, generally, if the quantum corrections are included, i?2 7^ 0. 

Now we will show that the quasiclassical approximation is exact in the soliton sector of sin-Gordon 
model. The structure of the perturbation theory is readily seen in the 'normal-product' form: 



N 



where 



and 



Ml) = E / UWmh ■■ e- lU ^-^ 2 ^e* s °^\r(q;u N )\ 2 :, 

Oun 



On 



J = J? • 



dr] 



Jv ■ 



Bun 



jx = / dt'Q{t - t')X(t') 



(K.41) 

(K.42) 
(K.43) 



with the 2iV-dimensional vector X = (£,??). In cq.( K.42j ) Q is the ordinary symplectic matrix. 

The colons in (K.41) mean that the operator j should stay to the left of all functions in the pertur- 
bation theory expansion over it. The structure (K.42) shows that each order over jx ( is proportional 
at least to the first ord er der ivative of un over the vari able conjugate to Xi. 

The expansion of (|K.4J|) over jx can be written [|105| in the form (omitting the quasiclassical 
approximation) : 

r N 2N r) 

M<t) = E/ IIw^mE dx^ PxiiUN)} > (K - 44) 



JV 
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where Px t {un) is the infinite sum of 'time-ordered' polynomials (see 1 105 ) over un and its derivatives. 
The explicit form of Px^un) is complicated since the interaction potential is nonpolynomial. But it 
is enough to know, see QK.42D , that 

P Xz (u N )~n tJ ^. (K.45) 



dX oj 



Therefore, 



R2(q) = 



(K.46) 



since (i) each term in (K.44) is the total derivative, (ii) we have (K.45) and (iii) un belongs to Schwarz 
space. 
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Figure 2: Keldysh time contour. 
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Figure 3: Niemi-Semenoff time contour. 
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Figure 4: Matsubara time contour. 



76 



10 



in 

N 



1 



10 



pp -> 

Vs 
Vs 



o 



1 GeV 
30.4 GeV 

52.6 GeV 



+ 



1 o 



o 



10' 



10 



® 



10 



-1 



10 



-2 



® 



10 



-3 







0.2 0.4 0.6 0. 



1.4 1.6 

z 



